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Abstract 


Modelling of blast furnace ironmaking has been carried out in four stages in the 
present work. Modelling of the raceway region enables to determine required 
blast humidity to maintain certain adiabatic flame temperature. Blast tem- 
perature, oxygen enrichment of the blast, pulverised coal injection rate, lime to 
coal weight ratio, ash in coke and degree of gasification of silica etc. have been 
considered as operating variables. Second stage of modelling includes heat and 
mass balances in the lower part of the blast furnace comprising of the hearth 
and bosh regions. It is aimed at finding the fuel rate to meet the heat demand 
of the furnace. Operating variables include the parameters mentioned for the 
first stage of modelling, shaft efficiency, heat losses and reserve zone temper- 
ature of gas and solid phase etc. Third stage of modelling is carried out to 
determine the amount, composition and temperature of outgoing gases. This 
helps to test the validity of the model by comparing the computed results with 
the plant data. Fourth stage of modelling allows to determine the amount and 
composition of intermediate slags formed in the bosh region. Injection of coal 
and lime would result in lowering of blast humidity to maintain the constant 
flame temperature. This, in turn, may decrease the fuel rate marginally. Re- 
placement ratio is found to be around 1.1 kg of coke per kg of coal injected. 
Injection of lime through tuyeres decreases the lime to silica weight ratio of the 
bosh slag. This, in turn, would help to improve the permeability of the bed and 
increase the productivity of the furnace. By maintaining lime to coal weight 
ratio of 0.2 at pulverised coal injection level of around 130 kg/thm. Computed 
basicity of the bosh slag decreases from 1.5 to 1.15 approximately and melting 
point of the bosh slag is lowered from 1690°C to 1540°C. Factors that favour 
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complete combustion of coal particles as well as complete dissolution of lime 
particles into the slag in the vicinity of the raceway must be studied carefully 
to take maximum advantages of Co-Injection of coal and lime through tuyeres 
of a blast furnace. 
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Chapter 1 


INTRODUCTION 


For many years, production of hot metal and steel has been subjected to 
increasing costs in the markets all over the world. Awareness concerning pol- 
lution of the environment has also been growing. The coke making process 
being connected with the blast furnace operation requires high expenditures 
regarding personnel and plant investment and is, moreover, polluting the en- 
vironment. Furthermore, coking coal supply is limited and in some cases must 
be purchased only by foreign currency. The scarcity of low ash coking coals 
coupled with the competition in blast furnace ironmaking, have forced iron- 
makers to innovate techniques to lower coke rate and produce high quality 
metal. Pulverised coal injection (PCI) is a notable innovation invented in 
the world. The pulverised coal injection in blast furnace is one of the fastest 
expanding development in the field of ironmaking because of the following: 

1. Pulverised coal is cheaper than coke resulting in improved economics. 

2. The capital cost of PCI is less than the coke oven plant on the basis of 
an equivalent unit of blast furnace fuel. The cost of PCI plant includ- 
ing preparation plant may be Rs. 1,500/- to Rs.2,300/- per annual tonne 
compare to Rs.3,500/- to Rs.4,000/- for the new coke oven plant[l]. 

3. A wide range of coal types from lignite to anthracite have been injected 
successfully with the ash level varying from 3 to 18 %[lj. 
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4. Conservation of coking coals will sustain the steel industry for longer 
periods. 

5 The coal injection decreases the raceway flame temperature allowing in- 
crease in blast temperature and oxygen enrichment thus intensifying the 
blast furnace operation and maximising the productivity. 

6. The coal injection leads to smoother furnace operation with less hanging 
and slipping as compared to all coke operation and it gives quick response 
to change in furnace operation. 

Injection of pulverised coal through the tuyeres into a blast furnace is being 
practised on a global scale these days to improve the stability and reduce the 
cost of blast furnace Ironmaking. It replaces the coke from the burden to 
a certain extent, and it is very significant especially from the view point of 
economics and environmental standards. Further reduction in the cost of steel 
making may occur due to lowering of the silicon and sulphur content in hot 
metal in the blast furnace. One of the means for simultaneous reduction of 
silicon and sulphur contents in pig iron is the flux injection into blast furnace. 
Lime injection is favoured to achieve following benefits: 

• Insitu fluxing of silica in the coke. 

• To lower basicity of slag in the bosh region. 

• To improve fluidity of the bosh slag. 

• To reduce flooding related problems. 

• To improve blast furnace productivity due to higher driving rates. 

• To reduce silicon transfer to metal. 

• Better control of sulphur in hot metal. 

With injection of pulverised coal and lime or dolime, the heat locally generated 
from pulverised coal combustion is transferred to nearby injected lime. 
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Today more than 100 blast furnaces all over the world are equipped with coal 
injection facilities. Tata steel commissioned a coal injection unit in March. 1991 
in the F blast furnace. The system developed by Klockener stahl techmk(KST) 
Fig 1.1 [2] was adopted. It is capable of an injection rate up to 150 kg/thm. 

The present work was undertaken to study the effects of coal and lime injection 
on the performance of the blast furnace. Stage wise modelling of the furnace 
based on material and enthalpy balances has been carried out to study the 
effects of various parameters on the fuel rate and bosh slag characteristics 
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Flarno Front Eqpt 



Figure 1.1: Flow diagram of the KST coal injection system 





Chapter 2 


LITERATURE REVIEW 


2.1 Preamble 


Recently coal injection in blast furnaces has been adopted for decreasing the 
coke rate to the lowest possible extent commensurate with the furnace thermal 
balance and raceway conditions. A summary of plants where the technology 
has been in use is included m Table 2.1. A variety of coal injection systems have 
been developed to handle various types of coals ranging in size from pulverised 
to granular [20]. 

Most of the systems developed have performed well, achieving 0.8 to 0.9 tons 
of coke replacement by 1.0 tonne of coal. The extent of coal injection is limited 
by the plant and process engineering considerations. Therefore, it should be 
ensured at the planning stage itself that the plant engineering permits high 
injection rates. System designs equipped with flow regulation using either 
individual or uniform tuyere systems are available today which can inject coal 
up to 200kg/thm. It is. thus, possible from the plant engineering point of view 
to inject coal with a uniform distribution at high rates. Process engineering 
limits may be imposed due to incomplete coal combustion during injection. 
This combustion is influenced mainly by coal quality, oxygen to carbon atom 
ratio during combustion and the blast furnace operating conditions. The basics 



2 2 Coal injection system 


6 


and pilot plant studies as well as the operating experiences of the furnaces 
where coal injection technology has been adopted are reviewed in this chapter. 


2.2 Coal injection system 

A complete coal injection system is shown in Fig 1.1. It comprises of the 
following units: 

1. Coal receiving and storage unit. 

2. Drying and pulverising unit. 

3. Injection unit. 

The details of each of the above units have been described by Taneja[2] and 
are summarised below. 

2.2.1 Coal receiving and storage unit 

Coal for injection is transported from the yard to the ground hopper by 
dumpers. The coal is conveyed to the raw coal bin by a vertical flexo-well 
conveyor. 


2.2.2 Drying and pulverising unit 

The coal being supplied generally has a moisture content of up to 12 % against 
a desired value of less than 1 %. Therefore, drying and pulverising of the coal 
are carried out simultaneously in the coal grinding mill which is injected with 
hot gases (Fig 1.1). Hot gases are usually supplied by a hot gas generator which 
uses coke oven gas, natural gas or oil as the heating medium. 



2.2 Coal injection system 


7 


Table 2.1: Present status of coal injection 


Plant 

Country 

Production 

(t/day) 

Coal Injection 
Rate (kg/ton) 

Year 

Inland Steel Co. 

U.S.A. 

8655 

150 


Arbed-Belvel[3] 

Luxemburg 

2126 

75.8 

1980 

Acesita[4] 

Italy 

5000 

170 

1981 

Cleveland [5] 

Iron works 

Netherland 

2.4 

259 

1981 

Oita works [6] 
Nippon Steel 
Corporation 

Japan 

8500 

52 

1982 

Kobe Steel Ltd. 
(Kakogawa)[7] 

Japan 

8000 

200 

1983 

SSAB Lulea[8] 

Germany 

9000 

140 

1983 

Shoudu Iron and 
Steel Company[9] 

China 

8564 

137.1 

1984 

Wuhan Iron and 
Steel Company[10] 

China 


175-220 

1986-87 

Sidmar[ll] 

Belgium 

4250 

117 

1988 

Mizushima No.4[12] 

Japan 

10658 

70-90 

1989 

Sollac 

France 

8814 

160 

1990 

Kimitsu No.4[13] 

Japan 

9567 

132 

1990 

Hoogovens[14] 

Netherland 

8228 

175-212 

1991 

Kokura[15] 

Japan 

— 

150 

1991 
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Plant 

; 

Country 

Production 

(t/day) 

Coal Injection 
Rate (kg/ton) 

Year 

Thvssen[16] 

Germany 

8915 

173 

1991 

Gary B.F 4[17] 

U.S.A. 

- 

160 

1993 

Taranto Works[18] 

Italy 

5000-6000 

200 

1994 
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To save energy, the waste gases from the stoves of the blast furnaces are also 
used for the drying operation. The inert nature of the gases helps in maintain- 
ing the oxygen content at low levels. 

Temperature of the hot gases entering pulveriser varies from 150 to 300 °C. 
The outlet temperature from the mill is 90 °C, and the moisture content of 
pulverised coal is 1 % or less. 

Grinding of coal in the mill is carried out by stream of hot gases flowing up- 
wards through nozzles positioned around the rotating table. The hot gases also 
remove the moisture from the coal and help to keep the partially ground coal 
circulating in the grinding zone. A static classifier with a system of adjustable 
vanes to vary the pulverised coal fineness is mounted directly on the top of the 
mill. Oversize particles fall to the bottom of the cone in a similar manner to 
that in a cyclone. The mill is designed to avoid fires and to contain explosions. 

The coal is fed into the mill by a drag chain conveyor which is provided with 
a variable frequency drive to adjust the rate of coal feeding. Different plants 
use hammer mills, horizontal ball crushers and vertical ball grinding mills for 
grinding as shown in Table 2.2, depending on installation and operating costs. 


2.2.3 Injection unit 

The injection unit mainly consists of pressure vessels and lance. 


2. 2. 3.1 Pressure vessels 

The coal from the fine coal silo falls into a vibrating sieve through the silo 
output valve. The undersize of this vibrating sieve falls into the pressure 
vessel through the pressure vessel inlet valve. The vent valve on the pressure 
vessel remains open during the filling of the vessel. When the vessel is filled all 
the aforementioned valves ( silo outlet valve, inlet valve and vent valve ) are 
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Table 2.2: Different types of grinding mills used in different plants 


Plant 

Drying Medium 

Grinding Mill 

KST 

Hot stove waste gas 

M.P.S. mill 

ARMCO 

Natural gas/Hot stove waste gas 

Bowl mill 

K-BFI 

Stove waste gas/Natural gas 

Vertical ball mill 

KOBELCO 

Hot Stove waste gas 

Roller mill 

TARANTO 

Blast furnace gas/Cowper gas 

Vertical Roller mill 

SHOUDU 

- 

Ball mill 

HOOGOVENS 

Natural gas 

— 

ACESITA 

— 

Vertical Ball mill 
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closed and the vessel is pressurised to 90 % of the nominal pressure so that it 
becomes ready for injection. Since there are two pressure vessels which work 
alternatinglv, the pressure vessel which is ready for injection waits for the zero 
weight or minimum weight signal from the other vessel. As soon as this signal 
is received, hydraulically operated dozing valve at the outlet of the pressure 
vessel opens and the fines are transported pneumatically through pressurised 
fluidisation for injection. It follows naturally then that the outlet/dozing valve 
of the other pressure vessel closes and the procedure of filling, pressurising and 
waiting is repeated for this vessel, though in order to maintain continuity 
during dozing, the dozing valve of the second pressure vessel closes only after 
the dozing valve of the first pressure vessel has opened by 25 %. The dozing 
valve allows injection to be performed according to a preset rate. The coal 
flows through a single pipeline into the coal distributor located near the cast 
house of the furnace. From the distributor, individual lines are taken to the 
tuyeres. These lines have the same number of bends and the same length to 
ensure almost equal distribution to the tuyeres and the deviation from tuyere 
to tuyere is within 5 9c. 

2. 2. 3. 2 Lance 

To enhance the efficiency of combustion, location of the injection lance is 
extremely important as it controls the mixing of coal particles with the blast. 
When the lance is located upstream from the tuyeres, the combustion efficiency 
improves but when the lance is moved too far upstream, there is no further 
improvement in combustion but ash will be deposited on the inner surface of 
the blow pipe and tuyeres (Fig 2.1) this causes clogging of the blast passage. 
The optimum location is dependent on the type of coal, the blast temperature, 
particle size and shape of the lance[21]. 

The coal stream immediately gets diverted along the direction of the hot blast 
irrespective of the inclination of the lance. Diameter of the lance used is 
generally 12 -15 mm. In most cases, the lance enters the blow pipe at a shallow 
angle, which maintains a narrow plume of coal in the tuyeres. A central plume 




Figure 2.1: Effect of lance position on combustion efficiency and ash deposition 
in the blow pipe 
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will minimise mechanical erosion of the tuyeres. 

Problem of ash deposition in the blow pipe with high ash content of coal have 
to be tackled by optimising the lance location. The desirable properties of coal 
described by Bhanu are[22] are as follows: 

• Ash content : Less than 10 %. 

• Fusion temperature : It should be greater than 1400 °C. 

‘ • Volatile Matter : It should be greater than 30 %. 

• Size • 80 % coal must pass through 200 mesh size. 


2.3 Safety precautions 

The following precautions are necessary to ensure safety of the system : 

1. Pulverised coal storage reservoir should be inert(non- oxidising). 

2. Build-up any of static electricity should be avoided i.e. all the equipments 
should be grounded. 

3. Sparks should be avoided. 

4. The oxygen and CO level should be monitored. Oxygen content should 
not exceed 8-10 %. Otherwise emergency nitrogen purging should be 
started. 
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2.4 Factors affecting coal combustion 

2.4.1 Volatile matter 

Low rank coals with high volatile matter contents give better efficiency of 
combustion, this may be a parameter which ultimately decides the suitability 
of any coal for combustion[23j. 

Rapid devolatisation of coal results in a higher degree burnout in the blow 
pipe. This in turn can increase sharply the temperature of gases in the blow 
pipe. Coals with low volatile matter have high carbon to hydrogen ratios and 
they give a high replacement ratio. But, generally high or medium volatile coal 
is preferred for injection purposes since their combustion efficiency is reported 
to be high. 

Tadehara et al.[7] have studied combustion efficiency of various coals having 
different volatile matter as a function of axial distance from tuyere, the results 
are shown in Fig 2.2. According to them coals with volatile matter approxi- 
mately 40 % has higher combustion efficiency as compared with coals having 
20 % volatile matter. 

Tar produced during pyrolysis of coal decomposes to soot during combustion. 
Hunty et al.[24] has reported that in the temperature range 1600-2000 °C 
about one third of mass of volatiles is transformed to soot. Tar and soot are 
undesirable materials in the blast furnace. Soot is less reactive than residual 
char and complete combustion demands a high air/fuel ratio and a long res- 
idence time not available during coal injection. Soot can cause an increased 
fuel rate and many serious operating problems such as blockage of raceways, 
decreased burden permeability, undesirable temperature distribution with in 
the furnace, and hanging and roiling of the burden. Thus generation of tar 
during pyrolysis of certain high volatile matter coals is detrimental. 

Although coals with volatile matter contents of about 25-30 % are preferable, 
the amount of volatile matter in coal should not exceed 35 % as there is a 
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chance of explosion. Further storage and transportation of such coals would 
be difficult owing to the problem of spontaneous combustion, particularly m 
dry months. 

2.4.2 Particle size 

The particle size of coal is important for achieving the desired degree of com- 
bustion in the raceway zone. If the combustion is not proper, then the particles 
would get swept into the burden and interfere with the gas flow. The upper 
size limit for pneumatic conveyance of coal through a tuyere is governed by the 
maximum lance diameter which can be inserted through the side of a blow pipe 
without unduly obstructing the flow of the hot blast. Generally, the smaller 
the grainsize of coal, the greater is the degree of combustion, due to more 
exposed surface area. It is reported by Wakuri that, if the particle size of the 
pulverised coal is 100 n m or less, and the volatile content is 25 % or more, 
then its combustibility is usually good[6]. According to O’hanlon at least a 
5T ratio of pipe diameter to particle diameter is acceptable for injection [25]. 

2.4.3 Blast parameters 

The combustion efficiency increases with an increase in blast temperature. The 
effect is more pronounced in the temperature range of 1000-1100 °C and wears 
down above 1150 °C, Shyuki, 1985[21] and Wakuri, 1983[6]. 

Pandey[26] has reported that under blast furnace conditions the combustion 
efficiency of coal improves more effectively by increasing the oxygen content of 
the blast than by raising the blast temperature. 

According to Ueno[27] the combustion efficiency increases slightly with increase 
in blast humidity, due to reaction of char and steam. The addition of steam 
theoretically lowers the flame temperature, it must be accompanied by increase 
in blast temperature or O 2 enrichment. 
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2.4.4 Oxygen to Carbon ratio 

Carbon is oxidised to CO gas in the raceway. This requires overall carbon to 
oxygen atom ratio of 1.0 in the gas phase. A part of carbon that is combusted 
comes from the coal and rest is from the coke. According to Peters et al [16] 
the limit of coal injection is reached when blast oxygen to coal carbon ratio is 
2.1. It is equivalent to coal injection rate of about 150 kg/ton hot metal. 

2.5 Effect of coal injection on blast furnace 
operation 


Typical operating data for some of the blast furnaces where the coal injection 
technology has been practised, are included in Table ‘2.3. The experiences 
gained by the plants from running of different blast furnaces are summarised 
below in this section. The most likely parameter to be affected by coal injection 
would be permeability. Studies by Goncalves[4] have shown that the perme- 
ability of the No.l blast furnace at Acesita experienced a downward trend with 
coal injection. Burden permeability was found to decrease when coal injection 
was increased from 115 kg/ton to 165 kg/ton. 

Silicon content in pig iron decreased from 0.72 % to 0.45 % when coal injection 
rate was increased from 30 kg/t to 150 kg/t at Acesita works. 

At Arbed-belval interesting observations have been reported by Koster[3]. 
Upon coal injection, there was no significant decrease in overall permeabil- 
ity though a slight decrease was reported in the bosh region. There was also 
no improvement in the productivity. However, the adiabatic flame tempera- 
ture dropped by 1.5 to 2 °K per gram of coal injected. One encouraging result 
of coal injection was significant improvement in the quality of hot metal with 
respect to low sulphur content, a phenomenon which is not clearly understood 
and can be a subject of further research. 

Nippon steel has successfully adopted PCI technology developed by Armco, 
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Table 2.3: 


Plant 

Hearth dia. 

(m) 

Blast 

Temperature 

(°C) 

Adiabatic 

Flame 

Temperature 
Temperature (°I<) 

Replacement 

Ratio 

Arbed-Belvel[3] 

8 

1132 

“ 

0.81 

Cleveland 

Iron works [5] 
(pilot plant) 

0.5 

870 

2175 


Oita works 
Nippon Steel 
Corporation[6] 

14.2 

1290 

2300-2350 

1.025-1.2 

Kobe Steel Ltd. 
(Kakogawa) [7] 

13.2 

1250-1300 


0.88 

Taranto Works [1 8] 


1163 

2200-2300 

0.8-0.88 

Hoogovens[14] 

13.5 

1250 

2150-2350 

0.9-0.95 
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Babcock and Wilocox for its use m large blast furnaces. Wakuri et al-[6] have 
reported about the performance of no.l blast furnace at Oita works which 
is based on technology used at Amanda blast furnace at Armco's Ashland 
works. The analysis of high speed filips showed that most of the coal injected 
burned up within a millisecond, i.e. before it had passed into the raceway. The 
injection rate was increased in three steps. First an injection rate of 30 to 40 
kg/t was attained, then it was increased to 50 to 55 kg/ 1, and later 80 to 100 
kg/t. Each injection rate was maintained for two months or longer to stabilize 
the operation. The most notable changes brought about by introduction of 
Pulverised coal injection were: 

1. Increase in productivity. 

2. Constant blast temperature of 1300°C. 

3. Decrease of blast humidity from 30 to 10 g/Nm 3 . 

4. Increase of ore/coke ratio from 3.5 to 4.25. 

5. Hanging of blast furnace was eliminated. 

Catastrophic tuyere failures, furnace hanging and variations in replacement 
ratio were major problems encounrered at the National steel and Bruceton 
furnaces as reported by Ostrowski[28]. Tuyere failure at Weirton No. 3 and 
No. 4 blast furnace due to coal injection was reported. It was observed in 21 
weeks time, that no. of tuyeres damaged increased from 50 without coal to 169 
with coal. Many a time slag and metal oozed out of the furnace through the 
tuyere cooler and the blow pipe. The severity of problem was further increased 
by reducing the air volume which resulted in filling up of blow pipes. 

Ostrowski[28] also reported that a significant hanging occurred at Bruceton 
blast furnace due to Pulverised coal injection. By reducing blast pressure, 
hanging could be overcome. Ostrowski[28] reported further that injection rate, 
blast temperature, blast oxygen, blast moisture and the degree of coal powder 
burning infront of tuyeres govern the replacement ratio. 
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Eichinger et al.[29] has reported the effect of coal injection on the performance 
of blast furnace in Taranto (ILVA). Alongwith 20 % reduction in coke consump- 
tion, productivity increased from 57 to 62 thm/m 2 /day with introduction of 
Pulverised coal injection. Initially oxygen enrichment of the blast was main- 
tained at high level due to operation safety but to abate the hot metal cost, 
this was reduced. 

Brouwer [14] reported that the performance of blast furnace in Hoogovens 
greatly improved, as a result of coal injection in terms of economy and hot 
metal quality etc. Blast furnace production capacity increased by approx- 
imately 10 % compared to an all-coke practice due to the decreased blast 
consumption. Coal injection levels up to 160 to 180 kg/ton hot metal could 
be handled without special burden and coke quality requirements. For injec- 
tion levels of 200 kg/ton of hot metal and above, specific problems occurred 
in running the furnace. 

Bocong et al.[9] have repoted on the experiences of Shoudu iron arrd steel 
company that varieties and qualities such as fineness of coal particles, the 
operational technique of the furnace and injection system influence replacement 
ratio. A high quality fine coal with low ash content accelerated the combustion. 
Moreover, the replacement ratio is also dependent on theoretical combustion 
temperature and gas flow distribution. 



Chapter 3 


PLAN OF THE PRESENT 
WORK 


The main objective of the present work is to analyse the effect of pulverised coal 
injection along with powdered lime through the tuyeres into the blast furnace, 
on various parameters like coke rate, top gas composition and temperature, 
allowable blast humidity, slag rate and basicity etc. The study has been divided 
into following four stages: 

1. RACEWAY REGION 

2. LOWER PART OF FURNACE 

3. OVERALL MASS AND HEAT BALANCE 

4. BOSH REGION 

Modelling of raceway region was aimed at finding the required blast humidity 
to maintain the adiabatic flame temperature for specified parameters. The 
lower part of furnace is considered to be more critical from heat point of view 
than the upper part. Heat and mass balance in the lower part have been 
used to determine the required fuel rate for certain raw materials and blast 
parameters. Overall heat and mass balance of the furnace has been used to 
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determine the top gas volume, composition and temperature. Fourth stage 
is an extension of second stage which uses only material balance to find the 
amount and composition of slag formed in the bosh region of the furnace. 
Mathematical model and the computer programme to find the solution have 
been tested using data from literature as well as data collected from G blast 
furnace of Tata Iron and steel company, Jamshedpur. Practical importance of 
the study is highlighted. 



Chapter 4 


MODEL EQUATIONS 


4.1 Introduction 


Mathematical model equations are developed in this chapter to determine the 
adiabatic flame temperature, fuel rate, slag rate, burden weights, top gas com- 
position and temperature. The work has been carried out in different stages, 
starting from the raceway region to the blast furnace taken as a single unit 
for making material and enthalpy balances. At each stage of work, certain 
assumptions and simplifications of the process are made. 


4.2 Mathematical Modelling of Raceway: 

4.2.1 Introduction 

Metal and slag temperatures in the hearth are determined to a great extent 
by the maximum temperature attained by the gases at the edge of raceway in 
front of the tuyeres of a blast furnace. Many controls in the blast furnace are 
necessarily based on keeping the flame temperature in the narrow temperature 
range by adjusting the operating parameters such as blast temperature, blast 
humidity and oxygen enrichment of the blast. A mathematical model, to deter- 
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mine the effect of operating variables on the calculations of flame temperature 
is described m this section. 

Raceway is a pear shaped high temperature zone m the vicinity of tuyeres. 
Preheated coke particles descend through the bosh to reach raceway. Pul- 
verised coal may be injected through the tuyeres along with the blast. Carbon 
in the coal and coke takes part in the combustion reaction with stoichiometric 
amount of oxygen coming with the preheated air blast. Raceway region may 
consist of an oxidising zone where carbon is first converted to carbon diox- 
ide and a reducing zone where carbon dioxide, thus formed as well as water 
vapours present in the blast in the gases react with carbon to form carbon 
monooxide and hydrogen. Any hydrocarbons present in the pulverised coal 
and undecomposed carbonates present in lime that are injected through the 
tuyere would dissociate in the raceway. Ash in coke and coal would combine 
with lime to form slag. However, a part of silica in ash may undergo a partial 
reduction at such a high temperature to form silicon monoxide gas. Ash may 
contain a little amount of haematite(Fe 2 03), which is reduced to iron. Rest of 
ash along with lime forms a part of slag. 


Species that enter and leave the raceway are shown schematically in Fig 4.1. 

4.2.2 Assumptions: 

1. There is no heat loss in the raceway. 

2. Incoming solid coke is preheated to 1500 °C before reaching the raceway. 

3. Fraction of silica being gasified is empirically assumed. 

4. Combustion efficiency of coal is 100%. 

5. Molten slag and hot metal descending through the bosh to the hearth, 
do not come in contact of the raceway. 



With \ 
Carrier Gas \ 


Preheated \ 
Air Blast 
with Moisture 


Preheated 
Solid Coke (T c ) 


Tuyere Gases (T^) 


CO H 2 N 2 SiO 



Tuyere Metal (Tf) 


Tuyere Slag (Tj) 


L. 


4.1. Schematic of Material Flow in the Raceway Region 
of the Blast Furnace . 
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6. Silica in the lime injected through the tuyeres is not reduced in the 
raceway. 

7 . Sensible heats of slag and molten iron are assumed to be independent of 
the composition of these phases and are determined by'empirical corre- 
lations. 

8. Sensible heat of the coal and the carrier gas are ignored. 

9. Heat of cracking of coal is determined from its ultimate analysis as fol- 
lows: 

• Any dissociable oxygen present in coal is assumed to be present 
with hydrogen as H 2 0. 

• Remaining amount of dissociable hydrogen is assumed to be present 
as CH 2 . 

• Enthalpy change due to transfer of sulphur and nitrogen for solid 
to gas phase are ignored. 

4.2.3 Overall Reactions 

Overall reactions occuring in the raceway may be described as follows: 

1. Oxidation of carbon by oxygen and moisture in the blast. 

C+(\)0 2 (g)^C0(g) (4.1) 

C 4- H 2 0(g) = CO(g) + H 2 (g) (4.2) 

2. Reduction of Si0 2 

S 1 O 2 + C CO(g) + SiO(g) (4.3) 

3. Calcination of undecomposed Lime 


CaCOz + C# CaO + 2 CO{g) 


(4-4) 
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4. Reduction of Iron oxide present m the coke and coal. 

Fe 2 0 3 4-3 C ^ 2 Fe + 3 CO{g) (4.5) 

Mathematical Modelling of Raceway is divided into two parts as follows: 

1. The blast parameters are specified and the model aims at finding the 
adiabatic flame temperature. 

2. The flame temperature reached by tuyere gases is specified and model 
aims at finding the required moisture in the blast for specified oxygen 
enrichment, blast temperature and level of pulverised coal injection rate 
etc. 

4.2.4 Model to determine adiabatic flame temperature 

One mole of dry blast consisting of oxygen and nitrogen is taken as the basis. 
All other parameters are assumed or determined per unit of dry blast. The 
blast includes any oxygen that is added to enrich the air. It is necessary to 
determine the amount of tuyere gases such as carbon monoxide, hydrogen, 
nitrogen and silicon monoxide as well as the metal and slag phases that are 
formed in the raceway, before making any heat balance to find the tempera- 
ture of the tuyere gases. The amount of coke that enters the raceway has to 
be determined as a function of operating parameters such as blast humidity, 
oxygen enrichment of blast, and amount of coal and lime that are injected with 
the blast. Model equations based on the balance of species such as nitrogen, 
hydrogen, carbon, oxygen etc. are described in Table 4.1. Symbols used are 
explained in nomenclature in the beginning of the thesis. 

There are 10 equations with 10 unknown variables. Variables n;v 2 , n# 2 , ncao, 
xiM g o, can be obtained directly from equation 1, 2, 7, 8 of Table 4.1. Equations 
3, 4, 5, 6 and 10 are solved simultaneously to determine n C o, n Sio , n Si o 2 , n Fe 
and W 2 . Then the value of n. 4 j 2 o 3 is obtained from equation 9 of Table 4.1. 
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Table 4.1: Equations based on Material Balance in the raceway region 


CO 

s: 

Balancing species 

Equation 

l. 

to 

+ (i - o«) + (in + wy/c, 

2. 

to 

n H 2 - ^ + X 

3. 

C 

I ( W 2 a+ wi Cl +0. 12 Wi 1 1 ) 

nco — 1200 

4. 

O 2 

_ 3(lViAvi-W2A2)a% W\0\ , 2{\V$Li) , op i T t r ^ n 

n CO — 1.6xl0 6 ^ 1600 10000 + + X + TISIO 

5. 

Si 


6. 

SiO 

I Wi M+W^MJEU 

n SiO — 6.0x10 s • 

7. 

CaO 

nn n - mLx 

fl CaO 56oo 

8. 

MgO 

4000 

9. 

AI 2 OZ 

~ _ W!Ai+W 7 AiM , 

'M/ 2 O 3 — 1.02x10® 1 1.02x10“ 

10. 

Fe 

2(WiA 1 +W' 2 a 2 )a? 
n /e ~ 1.6x10® 
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Heat balance of the raceway region gives the following equation 


n co Hf° + n Sl0 (tff 0 - H? l ° 2 ) + O t H T 0 \ + (1 - O e )H\ 

l H 2 0 


-j-xHuin + ~ 2 ~ 2 X 0.25 x (t c 


100 


298) + x H? c 


1200 


xHf-° + n H2 H T H {_ + n co H T c f 0 + n N ,H% 2 + n Sl0 H^ 0 

•F 

+n Fe (10.6(7) - 1873) + 280 x 56) + 


TjFeoOz 

ripe-ti f 


+(60 n Sl0 , + 102 71^,03 + 56n Ca0 + 40nA /sO )(0-25 x (7) - 1873) + 480) 

W; 3LlHf aC 


+ W iH Crack + 


10000 


(4.6) 


This is a nonlinear equation and it can be solved by the method iteration using 
a computer programme 

4.2.5 Model to determine blast humidity 


This model is also based on unit mole of dry blast. The model equations based 
on heat and mass balance are similar to those as described in section 4.2.4 and 
are summarised in Table 4.2. 

There are 8 linear equations with 8 unknown variables. Variable n^ can be 
obtained directly from Equation 1 of Table 4.2. Linear equations 3 to 8 of 
Table 4.2 are solved simultaneously to determine the unknowns x, W 2 , n co, 
risiCh W Fe and W s i ag . Then the value n # 2 may be obtained from equation 2 
of the Table. 


4.3 Modelling of the Lower zone : 

4.3.1 Introduction 

In this section, lower part of the blast furnace consisting of the bosh and 
hearth regions is considered. The lower part is separated from the upper 
part by a reserve zone where the temperature of gas is close to that of solid. 
Haematite(Fe 2 0 3 ) and magnetite (Fe 3 0 4 ) are assumed to be pre reduced to 
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Table 4 2: Equations based on Material and Enthalpy balances in the raceway 
region 


S.N 

Balancing 

species 

Equation 

1. 

iV 2 

n N2 = ™ + (1 - O e ) + Wi + 

2. 

H 2 

n H 2 - + X 

3. 

C 

(W 2 C2+W l Ci+0.l2W 3 L l ) 
n CO - 1200 

4. 

o 2 

~Z _ 3 (VK 1 4 1 +VV J 2 A 2 )a| W,0, , 2(W 3 Li) , no , „ , ~Z 

n CO — 1.6xl0 6 ^ 1600 ^ 10000 + + x + T15ZO 

5. 

Si 

~ (iViAi+lVaAaja?/ 

^zO “ 6.0x105 

6. 

Fe 

iVFe ~ 1.6x10® 

7. 

Slag 

■ ■ 

8. 

Enthalpy 

ncoHf° + n St0 (Hp° - Hf u >) + 0 e H l o \ + (1 - O e )H& 

|r r rT„ . (W 2 A 2 )(0 25)(Tb-298) . (W a C 3 )ffg 

w 100 l 1200 

= ®Sf»° + n„ 2 H% + ncolll! + n N ,H T N [ + n Si0 H T J 2 

+tV F .( (106 ^- le73) + 280) + 

w'oLi 7T Cac 

+t*W480 + (0.25) (Tf - 1873)) + W^Crack + Vodoo' 
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wustite(Feo.950) phase by carbon monoxide and hydrogen gases in the upper 
part. Wustite is partly reduced in the reserve zone by the ascending carbon 
monoxide and hydrogen gases. The amount of iron produced in the reserve 
zone would depend upon the following parameters: 

1 Reserve zone gas remperarure 

2. Reserve zone solid remperarure. 

3. Shaft efficiency. 


4.3.2 Thermodynamic consideration 


The indirect reduction reactions occurmg m the reserve zone may be written 
as follows: 


F & o 95O + 00{g) 0.95-F'e 4- C02(g) (4. 1 1 

Fe 0 95 O 4- 0.95 Fe + HiO( K g') (4.81 


The amount of metallic iron thar enters the lower part of the furnace is thus, 
related to the amount of reducing gases such as carbon monoxide and hydrogen 
that leave the bosh region. From the thermodynamic analysis we may derive 
the following: 


n F e = n 


ft CO 

1 -i- K± 7 


nh 2 

^ 1 + Ki 8 


(4.9) 


Where. 


4300 


logK i7 ( 4 606 (T Jt - 50) L14 ^ 


(4.10) 


_ __ f 43o0 \ 

logK t s = -( 4 606 ( Tji _ 50 ) + ■ ) 


(4.11) 


Reduction of the remaining amount of wustite would occur at high tempera- 
tures in the bosh. The reaction mechanism includes the gaseous reduction of 
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wustite and the carbon gasification reaction i.e.boudouard reaction which is 
described below: 


C0 2 {g) + C ^ 2 CO{g) (4.12) 

The direct reduction of wustite in the bosh may be described as follows. 


Fe 095 O + C ^0.95Fe + CO{g) (4.13) 

Besides direct reduction of iron ore in the bosh and combustion reactions 
occuring in the raceway as described in section 4.2. other changes oecuring in 
the lower part of the furnace may be described as below: 

1. Carburising of iron phase. 

2. Melting of metal phase. 

3. Melting of slag phase. 

4. Reduction of oxides of nonferrous elements such as silicon, manganese 
and phosphorus etc. 

5. Distribution of sulphur between slag and metal phases. 

The lower pan is very critical for smooth running of the furnace from heat 
and chemical point of view. Heat losses through the walls and bottom may be 
around 80 % of total heat losses from the furnace. Mass and enthalpy balances 
are used to determine the required fuel rate of the furnace, as a function of 
the operating parameters such as blast temperature, blast humidity, oxygen 
enrichment of blast, levels of coal and lime injection through tuyeres etc, as 
well as ass um ed parameters such as heat losses, shaft efficiency, temperatures 
reached by gases and solids in the reserve zone etc. Schematic to show incoming 
and outgoing species in the lower part is shown in Fig 4.2. 
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4.3.3 Heat and Mass balances 


To make heat and mass balance in the lower part, it is important to know the 
gangue to iron weight ratio in the iron bearing materials m the burden and 
the composition of the gangue. If iron ore and sinter percentage of the burden 
and their composition are known, one may get the gangue to iron weight ratio 
as follows - 


9r = 


IEj 9l 


(4 14) 


Where, 


e t is the percentage of species i in the burden, 
gj the weight fraction of gangue in species i. 
ft is the weight fraction of iron in species i. 

Gangue is assumed to consist of silica, alumina, mangnesia and calcium oxides. 
The following assumptions are made regarding the composition of hot metal: 


1. All P 2 O 5 in burden is reduced completely in the bosh to give phosphorous 
in metal. 

2. Partition coefficient of sulphur Ks is fixed using some empirical correla- 
tions. It is defined as follows: 

Ks = in slag 

9c S in Hot metal 

Weight fraction of sulphur in slag = Ks * fs 
Where, 

fs is the weight fraction of sulphur in hot metal. 

3. Recovery of manganese, R™ in hot metal may be known by intuition or 
otherwise. Manganese content of the burden species may not be known 
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Table 4.3: Equations based on Material and Enthalpy balances in the lower 
part of blast furnace 


S.N. 

Balancing 

species 

Equation 

1. 

n 2 

ft = ^ + (1 - 0.) + (IU + H 3 )Rc 

2. 

H 2 

1_ , 

P2 — 200 + X 

3. 

o 2 

PsSi , ncD _ 1.5(WiAi + W 26 A 2 )a? , , W 3 L 1 , n 

6900x2 U - D - r 7 — 1.6x10 s + 3200 ~ t ~ 10000 + 

+0.5x + 0.5P 4 + Pq( jg + ^ + =^ z ) 

4. 

Fe 

(0.95)(56)F s 5i , p / -I £ f £ £ \ 

6900 ~riT6^x ,/s Jmn JpJ 

= + 56 (P 8 + 0.95P 4 ) 

5. 

C 

f J P | X°jP (^iCi-f^C 2 + 0 . 12 U 3 Li) 

6. 

Slag 

n /1 Si r tv“ \ 1 0 /60 x 71 w /mn[l— H mn ) 

P5(.i ioo JaX-t) + ■rbK'toJsi 55 x iw 

= (^+^)(i- o-oi°D + 56g r (0.95P 4 + P 8 ) 

7. 

Shaft 

0.95^7^4.7 t 0.95p2^4.8 E> 

l-fi^.7 1 1+Jf4.8 “ 

8. 

Enthalpy 

PjHf + O e P^ 2 + (1 - O e )H£ + xH&o 

, (W 2 6A 2 )(0.25)(r 3e -298) . {W 2b C 2 )H^ c 
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= -xEf° + P 2 H T £ + P 7 Hjr + PiPjf 
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in most of the cases. MnO content of the slag may be determined by 
making manganese balance if manganese recovery as defined below is 
known. 


Rmn — 


(W mn ) HM 


(Wmn)j 


/ mnjslag 


(4.16) 


{WmnjHM = fmn x WrU 


(4.17) 


(1 = 5 x x Smn W HM (4.18) 

00 rlmn 

Burden calculations : 



Model equations based on the heat and mass balance in the lower part of the 
furnace by taking one mole of dry blast as the basis are summarised in Table 
4.3. There are 8 linear equations with 8 unknown variables, unknown variables 
are W 2 &, Pi, P2, P4, P5, P6, P7, Ps- Equations 3 to 8 of Table 4.3 are solved 
simultaneously to determine the variables W2 b and P 4 to Pg- Pi and P2 are 
determined directly from equations 1 and 2 of Table 4.3. Model equations to 
determine specific consumption rates of materials per 1000 kg of hot metal are 
summarised in Table 4.4. 
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Table 4.4: Calculation of different species per tonne hot metal. 


S.N. 

Species 

Unit 

Equation 

1. 

Coal 

Kg 

„ _ W, xlO a 

7i ~ pT~ 

2. 

Coke 


„ _ W?.bx.l0 6 

72 - PT~ 

3. 

Lime 


O, = Waxioj 

Pe 

4. 

Blast 


n — 103 

74 - 

5. 



Z _ lOixH 

^ P 6 

6. 

Moisture 
in Blast 

Kmoles 

„ __ lO^xx 

76 - p 6 

7. 

Iron 

Kmoles 

~ 10 3 x(0.95P 4 +P 8 ) 

qr _ 

8. 

Burden 

Kg 

o 

1 — ( 

II 

§ 

9. 

Iron ore 

Kg 

„ - HPxWor 

7°r — p 6 

10. 

Sinter 

Kg 

- _ 10 3 xW,„ 

— p 6 

11. 

Fe2<93 
in coal 


qixA\xa% 

710 = • T 6"1Q^ 

12. 

Fe203 
in coke 

Kmoles 

“ _ 92XAiXa? 

711 — l 6x10® 
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4.4 Overall Blast Furnace Modelling: 

4.4.1 Introduction 

In this section overall blast furnace has been taken as a single entity for con- 
sideration. Incoming and outgoing species are shown in Fig 4 3. The model 
is aimed at finding the temperature and composition of flue gases, that leaves 
the furnace for comparison and evaluation with actual data. Flue gases con- 
sist of carbon monoxide, carbondioxide, hydrogen, water vapour and nitrogen. 
Model equations are derived by taking one tonne of hot metal as the base. 
The blast rate, coal injection rate, lime injection rate and amount of burden 
materials such as coke, ore and sinter have been determined m section 4.3. 


4.4.2 Model equations 


Oxygen to Iron atom ratio in each burden species may be determined if total 
iron content as well as FeO content of the species is known. 


, hi (56 * 0.95) N x 56 

Fe 1 ' X ^ A x 69 )+ /,x6 9 


(4.19) 


Where, 

h 2 is the weight fraction of FeO in species i. 
Moles of Fe 2 0 3 being charged 


h = £(/* 


, 56 x 0.95.. 160 

FF )) x 


69 


112 


Moles of FeO being charged 




(4.20) 


(4.21) 


Overall oxygen to iron ratio in the burden may be calculated as follows: 




Moisture (T b ) Hot Metal (1873 °K ) 


Fig. 4.3. Schematic of Overall Material Flow in 
the Blast Furnace . 
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0 , US-J x 




56 


^ 56 


Carbon balance gives the following: 


(4.22) 


QiCj. . Q 2 C 2 0 12$ 3 I 


100 100 ' 100 
Oxygen balance gives the following : 

9iOi 3(?i-4i — 


"* 77174 ~~ 12(iVco + Nco 2 ) + 1000 x fc (4.23) 


1600 


+ 


1 6 x 10 6 


+ 96 + q&q 7 + 2 g 4 O e 


_1000 x (^1 + ^L + ^) 

V 00 1 err 1 / 


, g3^i x 0 02 
~~ 100 


28 55 62 

= iVco + 2Nco 2 + iVfl- 2 o 


(4.24) 


Hydrogen balance gives the following: 


Q1H1 , r , r , , . _ , 

<76 w + N H2 o (4.2o ) 

Relative amount of hydrogen, water vapours, carbon monoxide and carbon 
dioxide in the gas are influenced by the water gas shift reaction. 

H 2 0{g ) + CO(g) ^ C0 2 (g) + H 2 {g) (4.26) 


JV 

7* co 


C0 2 7 Y v N H 2 0 

z = A34 X 


N> 


H 2 


(4.27) 


The reaction being exothermic in nature is favoured at low temperature. How- 
ever the rate of reaction decreases exponentially with fall in temperature. For 
all purposes reaction may be ignored at temperature of less than 700K. 

There are 3 linear and 1 non linear equation. Which are solved by using a 
computer programme to determine Ntf 2 , N# 2 o> N^o and Nco 2 - 
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Nitrogen balance gives the amount of nitrogen in flue gas. 

QlNl +qi(l-O e )+R cg (q 1 +q 3 ) 


N, 


n 2 


2800 


(4.28) 


Overall heat balance is used to determine the top gas temperature. 
N co Hf° + N C0 2 Hf ° 2 + q,O t H T Q \ + 94 ( 1 - O e )Hj} 2 + q,H T H \ c 


NcoH'cS + N C o 2 Hqq 2 + N H2 H t £ + N N2 H% + N h , 0 H' 


1 9 


q z LiH? ac 

+qxH cr + q 4 H l0 + — ~ n J + q 5 h ls + 1000 H hm + q 2 M 2 H%° 


- 1000 ( 


1000 

f st Hf l °* fmnHf n ° f p H? 2 °s 


-( 


28 

1000 h 


)- (?6 - N H2 o)H^° 
+ <711 + gio)i?f 203 - i0 


(4.29) 


From heat balance, a cubic equation with one unknown T tg is obtained, which 
is solved by iteration. 


4.5 Fourth Stage of Modelling: 


Fluxes are either incorporated m the sinter or added directly to the furnace to 
form a slag of desired composition and melting range. Liquid metal and slag 
phases are separated out by gravity in the hearth or in the cast house. The 
composition of slag formed in the bosh may be diff erent from that of the cast 
slag because bulk of the ash in coke is released in the raceway region. 

Good permeability of the coke bed in the bosh region is very important for 
smooth running of the furnace. Presence of viscous slag in the gri^kiay affect 
the flow of ascending gases through the voids in the coke particles, and this 
in turn affects adversely the blast furnace productivity. Model equations are 
derived in this section to find the composition of slag formed in the upper part 
and lower part of the bosh. These slags have been called stack slag and bosh 
slag respectively in the present work for simplicity. 
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Table 4 5: Equations to determine the amounts and composition of slag in 
different parts of the furnace 


S.N 

Balancing species 

Equation 


Cast slag 

' 

1 . 

Si0 2 

(<?l Ai +<7 0 Ao )af 

fs _ q° T x°ri+q sn xsni + fh'n - 1 /„x6xl0 4 

100 28 

2 . 

CaO 

ns qorxor 2 -f-qsnXsri 2 

°2 “ 100 

3. 

MgO 

ns _ q~orxorz+q 3n xsnz 

°3 - 100 

4. 

> 

to 

O 

CO 

ns __ (qorXor4+q a nXsn4+(qiAi+q2A2)x£L) 

— 100 

5. 

MnO 

Cf = x x 1 - 29 x 100 

6 . 

Total Amount 

«“ - < c >' + °s + c i + c ‘* + x 

7. 

Basicity 

9 l 3= _ ^l 


Raceway slag 


8 . 

Si0 2 

1 ft = Mi + x £ + 

9. 

CaO 

D Q3Xlr3 

~~ 100 

10 . 

MgO 

E> q3xL4 

K * - :i ioo _ 

11 . 

AI 2 O 3 

ft = feft + X ^ 

12 . 

Total Amount 

514 = (i?i + i?2 + + -^ 4 ) 


Bosh Slag 


13. 

Si0 2 

Si = Ci -Ri + 

14. 

CaO 

B 2 ~ C 2 ^2 

15. 

MgO 

£ 3 = C| - R 3 

16. 

A1 2 0 3 

£4 = Q - £4 
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4.5 Fourth Stage of Modelling: 


44 


Stack slag may comprise of constituents present as gangue in the ore, sinter 
and raw fluxes. Coke particles remain more or less intact. In the bosh region, a 
part of carbon in coke is consumed via the boudouard reaction. Ash associated 
with carbon in the coke is added to the stack slag to form the bosh slag. 

Transfer of silicon and manganese to metal are assumed to be occuring in the 
vicinity of the raceway in the hearth region. Composition of the bosh slag, 
thus, would not be much influenced by the composition of the hot metal. If 
lime is injected through the tuyeres, corresponding amount of fluxes in the 
burden must be removed so that final slag composition remains more or less 
unchanged Model equations to determine the amount and composition of slag 
in the stack, bosh and hearth region of the furnace are summarised in Table 
4.5 



Chapter 5 


MODEL RESULTS 


The model equations derived in the present work have been solved by devel- 
oping a computer programme. Results of computer calculations are reported 
in this chapter. 


5.1 Raceway Model 


Two models described in section 4.2.4 and 4.2.5 for the raceway region were 
tested against each other as follows: 

1. For fixed humidity of the blast, the adiabatic flame temperature was 
determined by using model equations in section 4.2.4. 

! 2 By assuming the adiabatic flame temperature same as the value com- 

puted above, the required humidity is recalculated by using model equa- 
tion in section 4.2.2 and compared with assumed humidity in step 1. Any 
discrepancy could be attributed to error, either in model equations or in 
the computer programme and was sorted out before proceeding further. 
In the present work, model equations and computer programme were 
used for different sets of blast temperature, blast humidity and oxygen 
enrichment of the blast. 
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Computed flame temperatures are plotted against blast temperature at fixed 
blast humidity and oxygen enrichment in Fig 5 . 1 . Computed results are 
plotted against blast humidity in Fig 5 2 for different levels of coal injection. 
Results depend upon whether presence of ash in coke is taken into account or 
not and are influenced by the assumed weight fraction of silica in ash that is 
gasified. The presence of ash in coke has been ignored by many investigators. 
The empirical relationships as reported by Wegmann[30] and by Chung et 
al.[31] are as follows: 


T _ (0-9341 Tb + 8208TT - ff(2402 - 1.2177 Tb) + 94 76) 
F {1 + W + 2H) 


(5.1) 


and 


T f = 1559 + 0.839Tb - 6.033W S + 82.87 x ~ - 0.031 x ( lgggjgf £ ) 2 ( 5 . 2 ) 

Vj 5 V B 


where , 

T p = Raceway adiabatic flame temperature °C 
Tb = Hot blast temperature °C 
W = Oxygen in blast (Nm 3 /Nm 3 of dry blast) 

H = Humidity in blast(Nm 3 /Nm 3 of dry blast) 

Rpc = Pulverised Coal injection ratio (T/W) 

Vp = Wind Rate(Nm 3 /min) 

0 2 =oxygen enrichment(Nm 3 /w) 

W b =moisture in Blast(g/Nm 3 ) 

Computed values by using these relationships are plotted in Fig 5.1 for com- 
parison with the results of the model calculations in the present work. 



Flame Temperature (°C) Fuel Rate (kg/thm) 



Fig. 5.1. Fuel Rate and Flame Temperature versus 
Blast Temperature . 

a. Present Model with Ash (f = 1.0) 

b. ,, ,, (f — 0.5) 

c. ,, .» >> (f = 0.0) 

d . ,, ,, without Ash 

e . Wegmann’s Model 

f. Chung etal. Model 




Flame Temperature (°C) l“ ue ^ R°* e (kg/thm) 


540 



2300 


2250 


Blast Temperature =1100°C 
Oxygen Enrichment = 2 °/o 
Shaft Efficiency = 85 %> 


2200 


2150 


2100 


2050 


44.64 


2000 


10 12 14 16 


Blast Moisture (g/Nrrr) 


Fig. 5.2. Fuel Rate and Flame Temperature versus 
Blast Moisture . 
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Table 5.T Comparison of computed flame temperature values 


Parameter 

Unit 

Case 1 

Case 2 

Blast Temperature 

°K 

1300 

1300 

Blast Moisture 

(g/Nm 3 ) 

7.21 

7.21 

Oxygen enrichment 

% 

2.3 

m 

Computed 

Flame Temperature 

°K 



Present Model 
(Without Ash) 

°K 

2563 

2165 

Wagmann’s[30] Model 

°K 

2653 

- 

Chung’s[31] Model 

°K 

2799 

2489 

Reported value 
Davenport[32] 

°K 

2590 

2200 
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Empirical correlations used by investigators, in general, give higher values of 
the flame temperature, than the value computed from the model equations 
based on heat and mass balance in the raceway region. 

Davenport[32] determined the flame temperature for two sets of data by mak- 
ing material and enthalpy balances. The data and results obtained by him are 
summarised in Table 5.1 and they compare well with calculated values using 
the present model. 

Effects of changing the operating parameters such as blast temperature, blast 
humidity, oxygen enrichment and coal injection rate etc. on the calculation of 
the adiabatic flame temperature based on the model equations are however, of 
same order as determined from the empirical correlations. Results of finding 
the coefficient of these parameters are summarised in Table 5.2. 


5.2 Lower zone mathematical model 

As explained before, the model aims at finding the fuel rate of the furnace to 
meet the heat demand of the furnace based on certain assumptions. For given 
flame temperature, moisture content of the blast is determined using the race- 
way model. Model equations in the lower zone are solved simultaneously using 
the computer programme to obtain coke rate, slag rate and burden rate etc. 
Computed moisture in the blast and fuel rate are plotted against blast tem- 
perature at fixed flame temperatures and for two levels of oxygen enrichment 
of the blast of 2 % and 4 % in Fig.5.3. Results of the model calculations to 
find the fuel rate as a function of blast temperature for fixed moisture content 
in the blast are plotted in Fig 5.1, for shaft efficiencies of 0.65 and 0.85. The 
flame temperature is not affected by the shaft efficiency. 

Computed effects of individual parameters such as blast temperature, blast 
humidity, oxygen enrichment, coal injection, lime injection, hot metal com- 
position etc. on the determination of the fuel rate are summarised in Table 
5.3. Computed results of fuel rate and flame temperature determination are 
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Table 5.2. Comparison to show effect of various parameters on Flame temper- 
ature calculations ■ 


Parameter 

Unit 

Present Model 

Wegmann 

[30] 

Jin Kyung 
Chung[31] 

Blast Temperature 

+100°C 

71.0 - 72.0 

76.3 

83.9 

Blast Moisture 

+lg/Nm 3 

-(5.4- 5.5) 

-(3.5 - 3.6) 

-(6.0 - 6.05) 

Oxygen enrichment 

+1% 

43.0 - 44.0 

47.5 - 49.0 

63.5 - 67.2 

Coal injection 

lg/Nm 3 

-(1.84 - 1.9) 

— 

-(5.5) 

Ash in Coke 

Degree of silica 
gasification 
(f=0.0). 

(f=0.5) 

(6=1.0) 

+1% 

-(1.8 - 1.86) 

-(6.65 - 6.75) 

-(11.7-11.85 ) 

- 

- 

Lime injection 

LOI = 0.0 

LOI = 10.0 

lg/Nnr* 

-(1.15) 

-(1.49 - 1.51) 




CENTRA uBRMO 





Shaft Efficiency = 65 °/ 0 
Heat Loss = 2000 cat/moie 
Adiabatic Flame 
Temperature (°C) 


<0 

\*? 2 - 


Oxygen Enrichment 


1700 



cn 690 

-X 


O 670 


900 950 1000 1050 1100 1150 1200 

Blast Temperature (°C) 

3. Blast Moisture and Fuel Rate versus 
Blast Temperature . 
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Table 5.3: Effects of individual parameters on fuel rate Determination. 


' Parameter 

Unit 

Range 

Fuel Rate 
(kg/thm) 

Blast Humidity 

+1 g/Nm 3 

4-16 

(1.5 - 1.56 ) 

Blast Temperature 

+l°C 

900 - 1200 

-(0.18 - 0.22 ) 

Oxygen enrichment 

1% 

0.0 - 3.0 

-(3.5 - 4.5) 

Reserve zone Temperature 

+1°C 

850 - 1000 

(0.39 - 0.41) 

Heat Losses 

0.1(Gcal/thm) 

0.06 - 0.17 

28.6 - 35.9 

Shaft efficiency 

1% 

65 - 95 

-(2.8 - 3.5 ) 

Coal injection 

l(kg/thm) 

61.5 - 163.0 

0.282 - 0.284 

Hot metal [Si] 

0.1% 

' 

0.65 - 0.95 

2.51 - 2.53 

Hot metal [Mn] 

0.1% 

0.1 - 0.25 

0.32 - 0.38 

Hot metal [P] 

0.1% 

0.1 - 0.4 

1.92 - 1.96 

Ash in coke 

1% 

10- 19 

0.68 - 0.7 

Slag Rate 

+l(kg/thm) 

231 - 307 

-( 0.14 ) 

Lime Rate 

LOI = 0.0 

LOI = 10.0 

+l(kg/thm) 

0.0 - 26 

0.09 

0.18-0.19 
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plotted against lime to coal weight ratio at different values of LOI in Fig 5.4. 


5.3 Overall heat and mass balance 

Model equations were used to test the blast furnace data taken from the lit- 
erature.The operating data for Kimitsu, Taranto, SSAB Lulea and Tata blast 
furnaces are summarised in Table 5.4. 

Computed fuel rate and gas utilisation coefficient (rjco) are plotted against 
shaft efficiency at different heat losses for each blast furnace in Figures 5.5 to 
5.8. Gas utilisation coefficient ( 7 ico) have been determined from the computed 
top gas analysis as follows: 


%C0 2 

Vco ~ ( %CO + %C0 2 ) 


(5.3) 


In general fuel rate decreases with an increase in shaft efficiency and increases 
with an increase in heat losses from the furnace. Heat losses for the furnace may 
be 2 % to 5 % of heat demand, which is equivalent to 0.05 to 0.14 Gcal/thm 
or 1000-2000 Kcal/Kmole of blast. From the plots of the fuel rate versus shaft 
efficiency, one may determine the actual shaft efficiency for a particular plant 
corresponding to the reported fuel rate. This in turn allows determination of 
the CO utilisation in the shaft. Results are summarised in Table 5.5. 


5.4 Fourth Stage of Modelling 


It is an extension of modelling of second stage, it is aimed at finding the amount 
and composition ot slag formed in different regions of the blast furnace. e 
slag formed in the stack will consists mainly of the gangue present m the ore 
sinter and flux. Ash in the coke remains more or less intact and rt does not 
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Fig. 5.4. Flame Temperature 
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Table 5.4: Data of the different plants 


Parameter 

Units 

Kimitsu 

J apan[33] [34] [35] 

Taranto 
Italy [18] [36] 

SSAB Lulea 
Germanv[8l 

TATA 

India[l9] 

Blast parameters 



Blast Temperature 

°c 

1208 

1165 

1123 

1130 

Oxygen enrichment 

% 

1.28 

4.41 

0.0 

3.04 

Blast Humidity 

g/Nm 3 

41 

9.5 

16.4 

67.12 

Coke Data 


Coke Rate 

Kg/t 

372 

317 

464 

486 

Ash in coke 

(%) 

11.6 

9.77 

8.9 

18.7 

Fe 2 0 :i in coke ash 

HI 

5.1 

8.0 

■ 

mm 

Al 2 03 in coke ash 

(%) 

27.7 

28.0 

28.0 

30.0 

Si0 2 in coke ash 

(%) 

63.7 

64.0 

64.0 

60.0 

Moisture in coke ash 

(%) 

4.0 

4.0 

4.1 

4.0 

Coal/Fuel Data 


Coal Rate 

Kg/t 

125 

183.2 

0.0 

124 

Ash in coal 

(%) 

10.4 

9.12 

7.7 

16.48 

Fe 2 0,3 in coal ash 

mi 

8.0 

8.0 

8.0 

10.0 

A1 2 0 3 in coal ash 

(%) 

27.7 


28.0 

30.0 

Si0 2 in coal ash 

(%) 

63.7 

64.0 

64.0 

60.0 

Oxygen in coal 

m 

8.5 

5.2 

8.5 

5.3 

Nitrogen in coal 

(%) 

1.9 

2.08 

1.9 

1.62 

Hydrogen in coal 

(%) 

4.7 

3.23 

4.7 

3.9 

Carbon in coal 

(%) 

74.1 

80.4 

76.8 

71.8 
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Parameter 

Units 

Kimitsu 

Japan[33][34][35] 

Taranto 
Italy [18] [361 

SSAB Lulea 
Germany[8] 

TATA 

India[19] 

Sulphur in coal 

(%) 

0.5 

1.04 

0.5 

0.38 

Raw Material Data 

— — L 



Ore in Burden 

% 

50 

50 

50 

33.45 

Ore Analysis 




SiO -2 in ore 

(%) 

2.88 



1.13 

CaO in ore 





HH 

IVigd in ore 

m 

*“ " 



0.0 

ABO;* in ore 

(%) 

2.078 



1.49 

Pe in ore 


G3.75 



wm 

Sinter Analysis 


SiO— 2 in sinter 

m 

7.4 




CaO in sinter 

(%) 

14.1 




MgO in sinter 

(%) 

4.4 




AI 2 O 3 in sinter 


1.1 




FeO in sinter 

(%) 

7.7 




Pe in sinter 


mam 




Hot metal 
Composition 


C 

% 



4.67 

4.13 

Si 

% 



0.27 

0.74 

... 

Mn 




■■ 

mm 

P 

j mm 



mm 

U.2U2 


s 


0.051 


0.058 
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Fig. 5.8. TJ C0 and Fuel Rate versus Shaft Efficiency for 
G Blast Furnace of Tata Steel. 
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Table 5.5: Summary of the determination of Shaft Efficiencies 


Plant 

Fuel Rate 

Assumed Heat Loss 

Shaft efficiency 

Computed 


(kg/thm) 

(Cal/mole) 

(Gcal/thm) 

(%) 

Vco{%) 

Kimitsu[33] 

497 

1000 

0.05 

89.4 

55.8 



1500 

0.076 

92 

56.6 



2000 

0.11 

94 

57.5 

Taranto[18] 

500.2 

1000 

0.045 

77 

50.0 



1500 

0.069 

79.6 

51.0 



2000 

0.092 

81.8 

51.6 

SSAB 

464 

1000 

0.05 

77.8 

50.2 

Lulea[8] 


1500 

0.077 

'81 

51.4 



2000 

0.105 

83.8 

52.3 

TISCO[19] 

610 

1000 

0.056 

75.6 

47.0 



1500 

0.085 

78 

47.8 



2000 

0.115 

80.4 

48.7 
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contribute to the slag in the stack. In the bosh, a part of carbon in coke is 
consumed due to the direct reduction of iron ore as well as carburising of iron. 
Ash associated with the carbon in the coke might join the slag stream in the 
bosh. 

Cast slag comprises of the slag that descends from the bosh, ash that is released 
due to combustion of coal and coke in the raceway region and lime that is 
injected along with the coal. A part of silica that is released in the raceway 
may be gasified to give vapours of silicon monoxide in the tuyere gases. 


Si0 2 + C # SiO{g) + CO{g) (5.4) 

At lower temperatures in the vicinity of raceway, silicon monoxide vapours 
take part in the gas-metal and gas-slag reaction. 


SiO(g) + [C]^[Si] + CO(g) (5.5) 

SzO(g) + CO(g) ^ ( Si0 2 ) + C (5.6) 

Reduction of silica via slag metal reaction given below is considered to be 
sluggish and hence ignored. 


{Si0 2 ) + 2 [C] [Si] + 2CO(g) (5.7) 

* If silica could be reduced by carbon to give silicon in the bosh iron, lime to 
silica weight ratio in the bosh slag would be increased further. Computed 
amounts as well as lime to silica weight ratio in the cast slag, bosh slag and 
‘ stack slag are plotted against lime to coal weight ratio in Fig 5.9. Typical 
results of computation using 4 stages of model equations are included in Table 
5.6. 
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Table 5.6: Typical Results of computation using 4 stages of model equations. 



1 

2 

3 1 

4 

5 

Assumed Values 


RAFT (°K) 

2073 

2073 

2073 

2073 

2073 

Blast Temperature (°C) 

900 

900 

900 

900 

900 

Oxygen Enrichment (%) 

2.00 

2.00 

2.00 

2.00 

2.00 

Reserve zone 
Temperature (°K) 

1123.0 

1123.0 

1123.0 

1123.0 

1123.0 

Heat Losses (cal/mole) 

2000.0 

2000.0 

2000.0 


2000.0 

Shaft Efficiency (%) 

65.0 

65.0 

65.0 

65.0 

65.0 

Coal Rate (g/mole ) 

0.00 

2.0 

2.0 

2.0 

2.0 

Lime to Coal weight ratio 

0.00 

0.0 

0.1 

HI 

HI 

Computed Values 


Coal Rate (Kg/thm) 

0.00 

132.24 

132.22 

132.21 

132.20 

Blast Moist (g/Nm3) 

58.31 

24.12 

21.98 

20.91 

19.83 

Coke Rate (Kg/thm) 

692.35 

546.05 

544.39 

543.56 

542.73 

Replacement Ratio 

0.00 

1.11 

1.12 

1.13 

1.13 

Lime Rate Kg/thm 

0.00 

.00 

13.22 

19.83 

26.44 

Fuel Rate Kg/thm 

692.35 

678.29 

676.61 

675.77 

674.93 

Cast Slag Kg/thm 

316.39 

311.60 

311.32 

311.18 

311.03 

Stack Slag Kg/thm 

217.34 

217.74 

205.09 



Bosh Slag Kg/thm 

225.20 

226.58 

214.04 

207.77 

201.50 

Blast Rate Nm3/thm 

1592.34 

1525.57 

1521.40 

1519.31 

1517.22 

Heat Losses (Gcal/thm) 

0.14 

0.14 

0.14 

0.14 

0.14 
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5.4 Fourth Stage of Modelling 


Top Gas (Nm3/thm) 

2219.50 

2172.35 

2166.03 

f 

2162.86 

2159.69 

Top Gas Temp. (°C) 

321.56 

350.29 

352.85 

354.12 

355.40 

Heat Demand 
(Gcal/thm) 

2.7816 

2.6710 

2.6631 

2.6592 

2.6552 

CAST SLAG 
COMPOSITION 


Si02 (%) 

34.55 

34.06 

34.03 

34.01 

34.00 

CaO (%) 

33.17 

33.69 

33.73 

33.74 

33.76 

MgO (%) 

7.88 

8.01 

8.02 

8.02 

8.02 

A1203(%) 

22.72 

22.55 

22.54 

22.54 

22.53 

S(%) 

1.23 

1.23 

1.23 

1.23 

1.23 

FeO (%) 

0.28 

0.28 

0.28 

0.28 

0.28 

MnO (%) 

0.17 

0.17 

0.17 


■i 

Ca0/Si02 

0.96 

0.99 

0.99 

0.99 

0.99 

(CaO+MgO) /Si02 

1.19 

1.22 

1.23 

1.23 

1.23 

STACK SLAG 
COMPOSITION 


Si02 (%) 

23.92 

23.92 

25.26 

25.99 

26.76 

CaO (%) 

49.01 

49.01 

46.19 

44.64 

42.99 

MgO (%) 

11.65 

11.65 

12.22 

12.53 

12.86 

A1203(%) 

15.42 

15.42 

16.34 

16.85 

17.38 

Ca0/Si02 

2.05 

2.05 

1.83 

1.72 

1.61 

(Ca0+Mg0)/Si02 

2.54 

2.54 

2.31 

2.20 

2.09 
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BOSH SLAG 
COMPOSITION 


Si02 (%) 


31.58 

33.34 

34.30 

35.32 

CaO (%) 

47.30 

47.10 

44.26 

42.71 

41.06 

MgO (%) 

11.24 

11.19 

11.71 

11.98 

12.28 

A1203(%) 

18.90 

19.02 

20.13 

20.74 

21.38 

S (%) 

1.73 

H 

1.79 

1.84 

1.90 

FeO (%) 

.39 

.39 

.41 

.42 

.43 

MnO (%) 

.82 

.81 

.86 

.89 

.92 

Ca0/Si02 

1.51 

1.49 



m 

(Ca0+Mg0)/Si02 

1.87 

1.85 

1.68 

1.59 

1.51 

TOP GAS COMPOSITION 
(dry basis) 


CO (%) 

26.44 

25.48 

25.49 

25.49 

25.50 

C02(%) 

17.39 

17.27 

17.26 

17.26 

17.25 

Vco(%) 






H2 (%) 

4.44 

4.29 

4.13 

4.06 

3.98 

N2(%) 

51.73 

52.96 

53.11 

53.19 

53.27 































Chapter 6 


DISCUSSION 


6.1 Effect of blast parameters on performance 
of blast furnace 


Effect of individual operating parameters on the determination of blast humid- 
ity for maintaining the constant flame temperature and corresponding effects 
on the determination of fuel rate are summarised in Table 6.1. 

The effect of individual parameters such as blast temperature, blast humidity, 
oxygen enrichment and coal injection on the flame temperature and fuel rate 
calculation are summarised in Table 6.2. 


\ 

The fuel rate decreases with an increase in the flame temperature in all cases 
but the effect is the maximum when the flame temperature is increased by 
increasing the blast temperature and the effect is minimum when flame tem- 
perature is increased by increasing the oxygen enrichment of the blast.The limit 
to which the flame temperature could be increased might vary from furnace to 
furnace depending upon the design of the furnace and nature of raw materials 
used. Operating the furnace at flame temperature higher than the limiting 
value for a long time may lead to irregularities such as hanging, burning of 
tuyeres and hearth break out etc. 
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Table 6 1: Effect of various parameters on Fuel rate and Blast moisture deter- 
mination at constant Flame temperature 


Parameter 

Unit 

Range 

Blast 

Moisture 

(g/Nm 3 ) 

Fuel Rate 

(kg/thm) 

Blast Temperature 

+1°C 

900 - 1200 

0.145 - 0.166 

-(0.034 - 0.037) 

Oxygen enrichment 

1% 

0.0 - 3.0 

8.13 

7.5 - 9.0 

Reserve zone 
Temperature 

+1 °C 

850 - 1000 

0.0 

0.4 - 0.43 

Coal injection 

l(kg/thm) 

66.0 - 166.0 

-0.26 

-(0.106 - 0.107) 

Ash in coke 

1% 

10-19 



(f=0.0) 



0.22 - 0.25 

0.4 - 0.5 

(f— 0.5) 



1.0 - 1.05 

0.88 - 0.9 

(6=1-0) 



1.7- 1.75 

1.8 - 2.0 

Lime Rate 

+l(kg/thm) 

0.0 - 26 



LOI = 0.0 



- 0.14 

-(0.13 - 0.14) 

LOI = 10.0 


1 

- 0.190 

-(0.11 - 0.12) 
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Table 6.2: Effect of blast parameters on Flame temperature and Fuel rate 
determination 


Parameter 

Unit 

AT f 

AW fuel 

fuel 

±T, 

Blast Temperature 

+ 1°C 

0.7 

-0.21 

-0.31 

Blast Humidity 

+ 1 g/Nm 3 

-5.5 

1.5 

-0.27 

Oxygen Enrichment 

+ 1% 

43.5 

-4.0 

-0.1 

Coal Injection 

+ 1 g/Nm 3 

-1.9 

0.46 

-0.24 
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Result shows that fuel rate tends to decrease at the rate of around 0 035 kg 
per unit increase in blast temperature, if humidity is increased simultaneously 
at the rate of around 15 g/Nm 3 or (0.25 kg/ton) to maintain the same flame 
temperature. If blast humidity was maintained at same level and flame tem- 
perature was allowed to increase, fuel rate would decrease at the rate of 0 21 
kg per unit increase in blast temperature. If the flame temperature could be 
maintained at the constant level by decreasing the oxygen enrichment of blast, 
at the rate of 0.017 % per unit increase in blast temperature instead of increas- 
ing the oxygen enrichment of blast, fuel savings would be the order of 0.133 kg 
per unit increase in blast temperature. If coal is injected, steam injection must 
be lowered at the rate of 0.3 kg per kg of coal injected in order to maintain the 
constant flame temperature. Fuel savings are computed to be around 0.1 kg 
per kg of coal injected and it gives a replacement ratio of 1.1. If coal injection 
rate is accompanied by increase of 1.7 °C in blast temperature, to maintain 
same temperature instead of decreasing the humidity, fuel savings are around 
0.16 kg, it gives coke replacement ratio of 1.16. On the other hand, if coal 
injection is accompanied by increase in oxygen enrichment of 0.03 % of blast 
for maintaining the flame temperature, fuel rate increases by 0.16 kg, to give a 
coke replacement ratio of 0.84. These results tend to highlight the fact that the 
oxygen enrichment should be practised when the maximum limit of increasing 
the blast temperature and decreasing the blast humidity have been reached. 
If coal injection beyond the limiting conditions is not accompanied by oxvgen 
enrichment of blast, there will be fall in the flame temperature. It may lead to 
chilling of hearth, incomplete combustion of coal particles and non dissolution 
of lime particles which ,in turn, would make it difficult for smooth working 
of the furnace. Oxygen enrichment helps to reduce the amount of gases pro- 
duced in the raceway and bosh region and it will have beneficial effects on the 
furnace productivity. Specific consumption rates of the blast for two cases are 
compared below in Table 6.3. Higher proportion of oxygen in the blast would 
also favour complete combustion of coal particles in the raceway as mentioned 
by K. H. peter[16]. 
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Table 6.3: Specific consumption rates of the Blast with oxygen enrichment 
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6.2 Reserve zone temperature 


It has no effect on the flame temperature calculations. Reserve zone tempera- 
ture might be decreased >vith an increase in the reactivity of coke. Utilisation 
of CO gas for indirect reduction in the reserve zone increases with a decrease 
m the reserve zone temperature. Results of model calculations give fuel saving 
of 0.4 kg/thm per unit decrease in the reserve zone temperature. Hence to get 
optimum results, coke of high reactivity should be preferred, resulting m lower 
reserve zone temperature and better utilisation of gases. 


6.3 Shaft efficiency 

It has got no bearing on the flame temperature calculations. The shaft effi- 
ciency is related to attainment of equilibrium for the reduction of wustite by 
gases such as CO and H2. Fuel rate tends to decrease at the rate of as much 
as 3 to 4 kg per unit percentage increase in the value of shaft efficiency. In 
modern large blast furnaces using well prepared burden material, the shaft 
efficiency may reach 100 % and there may not be much scope for reduction of 
fuel rate. Presence of hydrogen content in the gases in the bosh region tends to 
increase the shaft efficiency as hydrogen molecules being smaller diffuse faster 
into the iron ore than carbon monoxide. Hydrogen content of the gas increases 
with a increase in steam injection rate and coal injection rate. All efforts are 
required to improve the reducibility of incoming iron ore and sinter and other 
iron bearing species in the burden for getting shaft efficiency close to 100 %■ 


6.4 Heat Losses 

Heat losses in Gcal/hr. are estimated empirically as follows: 


ffioss 1.9 d/i 
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0.5 Lime injection 


whom (I/, is the diameter of health 

hoi laige and modern blast lurnaces heat losses are estimated to be less than 
0 1 Gcal/t.hm Increase in heat losses bv 0 01 Gcal/thm increases the fuel rate 
by 3 kg/thm, ihe extent of heat losses may vary from furnace to furnace. 
Uncertainty regarding heat losses may pose problems in checking the validity 
of the model. 


6.5 Lime injection 


If lime, is injected without changes in blast parameters such as blast tempera- 
1, me, blast humidity, oxygen enrichment etc., calculations show that fuel rate 
must increases by 0.09 kg per kg of fullv calcined lime and flame temperature 
decreases at the rate of 1.1 °C per g/Nm 3 of lime in the blast. If lime contains 
undecomposed carbonates so that loss on ignition becomes 10 %, decrease in 
flame temperature is 1.5 °C per g/Nm 3 of lime and fuel rate must increase at 
the late of 0.18 kg per kg of lime. If lime injection is accompanied by a de- 
crease in humidity or increase in blast temperature so that flame temperature 
remains unchanged, there may be fuel saving of 0.13 kg per kg and 0.11 kg per 
kg of lime at 0.0 and 10.0 % LOI in lime respectively. Results are plotted in 
Fig 6 1. One may conclude that heat balance consideration in raceway and in 
the lower zone of the furnace are not much affected by injection of lime up to 
the limit of lime to coal weight ratio of 0.2. It has, however significant effect on 
the amount and composition of inter mediate slags formed in the bosh region. 
Results of finding the amount and composition of bosh and stack slag have 
been summarised in Table 5.6. One may estimate the melting point of slag by 
i (’foiling to Al/).s - C!aO - Si(U tommy diagiam[37] which is shown in Fig 6.2. 
Presence of magnesia in the slag is considered to be equivalent to presence of 
lime for estimating the melting point using the ternary diagram. The amount 
and composition of cast slag have been assumed to remain independent of lime 
to coal weight ratio in the injcctant. Estimated melting points of the bosh and 
stack slag for different lime and coal weight ratios are plotted in Fig 6.3. 



Fuel Rate (kg/thm) Blast Moisture (g/Nm ) 
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Weight °/o AI2O3 — 


Figure 6.2: Ternary Diagram of A1 2 0 3 - CaO - Si0 2 



Lime to Coal Weight Ratio 

Fig. 6.3. Melting Point and % Al 2 0 3 versus Lime to Coal 
Weight Ratio . 
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Injection of lime through the tuyeres has beneficial effect of lowering the melt- 
ing range of the bosh slag and stack slag. Presence of oxides of iron and 
manganese in the slag may lower the melting range of these slags further. If 
melting range of the bosh slag and its fluidity are improved by injection of lime, 
the furnace might accept increased wind rates without running into problems 
such as flooding etc. The furnace productivity would, thus be enhanced. Re- 
activity of lime is very important. It refers to the rate at which it becomes part 
of the liquid slag at any temperature. If lime of poor reactivity is injected, lime 
particles may not dissolve in the slag and ultimately find their way to the bosh. 
The permeability of the coke bed would be adversely affected, and capacity of 
the furnace to accept wind may be retarded at the cost of the productivity. 
Efforts should be made to operate the furnace at slightly higher flame temper- 
ature when lime and coal particles are co-injected through tuyeres. This will 
ensure complete combustion of coal particles as well as rapid dissolution of the 
lime. Such conditions favour reduction in coke rate, an increase in the pro- 
ductivity and overall lowering the cost of hot metal. Injection of lime into the 
raceway could favour the following in addition to heat balance considerations: 

1. Lime combines insitu with silica, that is released due to combustion 
of coal and coke particles in the raceway region. This in turn retards 
gasification of silica to silicon monoxide. As transfer of silicon to metal 
is attributed to formation of silicon monoxide, silicon content of metal 
may decrease with increase in injection of lime through the tuyeres. 

2. Injection of lime increases the sulphide capacity of the slag in the vicinity 
of the raceway region. This will reduce the extent of sulphur pickup by 
metal phase in the bosh region. 

Benefits of lime injection on lowering of silicon and sulphur level in hot metal 
have been confirmed by Komatsu et al[38]. Results of their findings are sum- 
marised in Table 6.4. Scientific and planned plant trials are needed to assess 
merits of simultaneous injection of coal and lime through the tuyeres of a blast 
furnace along with the blast under indian conditions. The study must be aimed 
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Table 6.4: Operational results with flux injection at Kokura No.2 BF[38] 


Items 

Base Period 

Dolomite 
= lOkg/THM 

Dolomite 
= 15kg/THM 

Productivity (T/D) 

3581 

3554 

3567 

Coke rate (kg/THM) 

464.1 

470.5 

467.7 

PC rate (kg/THM) 

23.6 

22.9 

22.8 

Dolomite rate (kg/THM) 

0.0 

9.0 

13.5 

Blast temperature (°C) 

1178 

1173 

1181 

Blast moisture (g/Nm 3 ) 

33.3 

| 

33.8 

34.3 

Flame temperature (°C) 

2303 

! 

2294 

2296 

HMT (°C) 

1489 

1491 

1492 

Si (%) 

0.40 

1 

0.38 

0.38 

S (%) 

0.020 

0.019 

, 

0.018 

MgO (%) 

7.49 

7.17 

— - - 

7.29 

CaO/Si0 2 (-) 

1.23 

1.26 ! 

: 

1.27 
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at finding the operating parameters that give optimum results for a particular 
furnace. 



Chapter 7 


CONCLUSIONS 


Mathematical modelling of the blast furnace has been carried out in stages 
in the present work. The raceway region has been modelled to predict the 
adiabatic flame temperature for given blast parameters. Another model is de- 
veloped to determine the required blast humidity to maintain certain adiabatic 
flame temperature. Presence of ash in coal and coke and gasification of silica 
m the raceway have been incorporated into the model. Important findings of 
the study are as follows: 

1. Flame temperature increases by 0.7 per unit increase in blast tempera- 
ture. 

2. Flame temperature decreases by 5.5 °C per g/Nm 1 2 3 4 5 of moisture in the 
blast. 

3. Flame temperature increases by around 43°C per unit % increase in 
oxygen enrichment of the blast. 

4. Flame temperature decreases by about 1.9 °C per 1 g/Nm 3 injection of 
coal. 

5. The flame temperature decreases by 1.8 to 11 °C for increase of 1 % ash 
in coke depending upon whether silica is gasified or not. 
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6. Injection of lime of lg/Nm 3 decreases the flame temperature by 1.1 to 
1.5 depending upon presence of undecomposed carbonates up to 10 9c. 

Heat balance in the lower zone of furnace has been used to determine the 
fuel requirement of the furnace. Important parameters which affect the heat 
balance include the shaft efficiency, heat losses, gangue to iron weight ratio 
in the burden, blast parameters, coal injection and coke reactivity etc. The 
model may be used to evaluate the reported data of a plant and the scope for 
improvement. Analysis of the data from the literature for a number of modem 
and large blast furnaces shows that the shaft efficiency could be around 80 to 
90 %. Study shows that efforts should be made to operate the blast furnace 
at as high flame temperature as possible without causing irregularities such as 
hearth breakout and tuyere burning. Oxygen enrichment of the blast should 
be practised when the maximum limit of increasing the blast temperature and 
decreasing the blast humidity have been reached. 

Amount, composition and temperature of outgoing flue gases are determined 
in the third stage of modelling by making overall mass and heat balance of 
the process. Access to reliable plant data regarding top gas composition and 
temperature are required for the validity of the model. 

The last stage of modelling has been aimed at finding the composition and 
amount of the intermediate slags formed in the upper and lower bosh parts 
of the furnace. The melting points of the slags are determined using ternary 
Al 2 0 3 -Ca0-Si0 2 diagram. Melting point of the bosh slag decreases from 1690 
to 1540°C, when lime is injected along with coal to have lime to coal weight 
ratio of 0.2 at coal injection rate of 2g/mole or 132 Kg/thm. 
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APPENDIX 


6a 


COMPUTER PROGRAMME 
FOR SOLVING 4 STAGE MODEL EQUATIONS 

REAL Wck, Wl,0e,02,N2,Hdem( 1 0) 

REAL Acl,Ack,Aclsi,Acife,Aclal,Acksi,Ackfe,Ackai 

REAL OcI,NcI,Ccl,Hcl,FCck,Lsi,I .ca,Lmg,Lal,Lcaco3,f 

REAL Mn2 

REAL Hcrack 

REAL Sh2,Sn2,Sco,Ssio 

REAL Tf,Tc,Tbosh 

REAL Wcl( I0),Wlcac 

REAL Tb,a( 1 0, 1 0),c( I ()),y( 1 0) 

REAL Mh2o( 1 0),or( 1 0),sr( 1 0) 

REAL cps,cpsio,cpco,cph2,cpco2,cpn2,cpo2tb,cpn2tb,cpc 
REAL ac I ,bc 1 ,cc I,ac2,bc2,cc2,ac3,bc3,cc3,ac4,bc4,cc4 

REAL ac5,bc5,cc5,ac6,bc6,cc6 

REAL ac7,bc7,cc7,ac8,bc8,cc8,ac9.bc9,cc9 
REAL Hffeo,Hfsi,Hfmn,Hfp,Hfh2o,Hfcac,Hslag,Hhm,cpfe 
REAL fc,fsi,finn,fp,Mbn2,Mbh2,Hfco,HfFo3,Hfsio 
REAL KI,K2,Mfe2 

REAL u( 1 0, 1 0),v( 1 0, 1 0),vr( 1 0),ur( 1 0),p( 1 0),bs( 1 0),s( ! 0) 
REALurl,ur2,dl,d2,H4 

REAL CL,CK,BL,SL,LM,Mst;r co, [ co2,TIi2,Th2o,Tn2,FEcI,FEck,FEot 

REAL FEin,CAIm,02ty,N2ty,Hin,H I ,H2,H3,Stpg,Hfco2 

REAL C 1 ,C 1 2,C22,C2,C3 ,C4,C5,C6,C7,C8 

REAL root,cf,ffes,Tg,Tgp,si,ca,mg,al,rmn,fs,ks 

REAL WcoaI( 1 0), Wcoke( 1 0),Blast( 1 0),WfueI( 1 0), Wslag( 1 0) 

REAL ph2o( 1 0),WIime( 1 0),qa( 1 0),qb( 1 0),qc( 1 0),qd( 1 0),qe( 1 0) 

REAL Tgt(I0),tt(10),t,rr(I0),flIs,Loy,xms,msck,gr0 
REAL por,h( IO),FOin,FEOin,H5 
REAL Wcs( 1 0), Wbs( 1 0),lmpt( 1 0),eta,ffeos 
REAL bca( 1 0),bsi( I0),bmg( 1 0),bal( ID) 

REAL slca( 1 0),slmg( 1 0),slsi( 1 0),slal( I0),slmn( 1 0) 

REAL bs I ( 1 0),bs2( 1 0),slfe( 1 0),sls( 1 0 ) 

REAL coal( 1 0),bs3( 1 0),bs4( 1 0),bs5( 1 0),bs6( 1 0) 

REAL Wsit,Wslib( 10),sit,cat,mgt,alt 

REAL sib( 1 0),cab( 1 0 ),mgb( E0) ? alb(I0),sb(I0),feb(10),mnb(I0) 

REiAJ, nn2,nh2,ncao,nmgo,nA12o3,siuib,mnuib 
REAL rft( 1 0),bt( 1 0),oxy( 1 0),Gtm( 1 0),PGtm( 1 0) 

REAL HIoss( 1 0),efcy( 1 0),rzt( 1 0),L4l,LI ,tpc,to,tl 1 


INTECifiR j,zt,i-3 

character yy I*20,yy2*20,yzl*20,yz2*22,yz3*20,yz4*20,yz5*20 



character yz6*20,yz7*20,yz8*20,yz9*20,yzl0*20,yzl 1*20 
chai actcr y/. 1 2*20,yz 1 3 *20,yz 1 4*20,yz 1 5*20,yz 1 0*20,yz 1 7*20 
character yz 1 8*20,yz 1 9*20,yz20*20,yz2 1 *20,yz22*20 
character yz23*20,yz24*20,yz25*20,yz26*20,yz27*20 
character yz28*20,yz29*23,yz30*24,yz3 1*24 
charactci yz32*22,yz33*24,yz34*24,yz35*23,yz36*23,yz37*23 
open(20,file~'datatb.in') 
open(2 1 ,fiIe“'result.out') 

read(20,*)AcI,AcIfe,Aclal,Aclsi,Ocl,NcI,Hcl,Ccl,ScI 

read(20,*)Ack,Ackfe,Ackal,Acksi,FCck,xms 

read(20,*)si,ca,mg,aI,T,oy,f 

iead(20,*)Ha,Hb,Hc,Ild,He,calv 

read(20, * )ac 1 ,bc I ,cc I 

read(20,*)ac2,bc2,cc2 

read(20, * )ac3 ,bc3 ,cc3 

read(20,*)ac4,bc4,cc4 

read(20,*)ac5,bc5,cc5 

read(20, * )ac6„bc6,cc6 

read(20, * )ac7,bc7,cc7 

read(20,*)ac8,bc8,cc8 

read(20, * )ac9,bc9,cc9 

read(20,*)Hih2o J Hffe,Hfsi,Hfmn,Hfp,Hfcac 

read(20, * )HfFeo,HsIag,Hhm,Htco,Hfch2 

read(20,*)Hffo3,Hfsio,Hfco2,Hvp 

read(20,*)fcp,fsip,hnnp,fpp,fsp,rmnp,ks 

read(20,*)rb,lf,Tbosh,Hls,jt,dwcl,zt 

read(20,*)por,eta,Oe 

read(20,*)oi ( 1 ),or(2),or(3),or(4),or(5),or(6),or(7) 
read(20,*)sr( I ),sr(2),sr(3),sr(4),sr(5),sr(6) 
read(20,*)lmpt( 1 0),ffeos 

write(2I,*)Ad,Aclfe,Aclal,Aclsi,OcI,Ncl,Hcl,CcI,Scl 

write(2 1 ,*)Ack, Ackfe,Ackal,Acksi,FCck,xms 

write(2 I ,*)si.ca,mg,al,Loy,f 

write(2 1 ,*)Ha,Hb,Hc,Hd ; He,calv 

write(2 1 ,*)ac I ,bc I ,cc I 

write(2 1 ,* )ac2,bc2,cc2 

write(2 1 ,*)ac3,bc3,cc3 

writc(21 ,*)ac4,bc4,cc4 

write(2I,*)ac5,bc5,cc5 

write(2 1 , * )ac6,bc6,cc6 

write(2 1 ,*)ac7,bc7,cc7 

write(2 1 ,*)ac8,bc8,cc8 

write(2I,*)ac9,bc9,cc9 

write(2 1 ,*)Hfh2o,HfFe,Hfsi,Hfmn,Hfp,Hfcac 

write(2I,*)Iinco,HsIag,Hhm,Hfco,Hfch2 



So 


write(2 1 , + )I IfTo3 ,Hfsio,I Ifco2,Hvp 
write(2 1 ,*)fcp,fsip,fmnp,fpp,fsp,rmnp,ks 
write(2 1 ,*)Tf,Tbosh,IIIs,jt,dwcl,zt 
write(21,*)por,eta,Oe 

write(2 1 ,* )or( I ),or(2),or( 3 ),or(4),or(5),or(6),or(7) 
wi ite(2 1 ,* )si ( I ),sr (2),st (3),sr(4),sr(5 ),sr(6) 
write(2 1 , * )lmpt( 1 0),fFeos 
yyl ' Mast Moist r/Niii) ' 
yy2~' Blast I emp C = ' 
yzl~' Coke Rate Kg/thm = ' 
yz2~'Coal Rate Kg/thm - ' 
yz.!- Fuel Rate Kg/thm = ’ 
yz4- Cast Slag Kg/thm = ' 
yz5- Blast Rate Nm3/thm= * 
yz6=' Lime Rate Kg/thm = ’ 
yz7— CO (%) 
yz8- C02(%) 
yz9=’ H2 (%) 
yzI0='N2(%) 

yzl 1=' Top Gas Temp (C) = ' 
yz 1 2- Top Gas Nm3/thm = f 
yzl 3- Replacement Ratio = ' 
yzl4- Eta-CO (%) 
yz 1 5- New Gangue Ratio = 1 
yzl6--'Si02(%) 
yz 1 7~' CaO (%) ~ * 

yz 1 8- MgO (%) 
yzl 9- A1203(%) 
yz20- FeO (%) 
yz2I- MnO (%) 
yz22- Ca0/Si02 
yz23=' (CaO t Mg0)/Si02 = * 

yz24= r S (%) 
yz25“' Stack Slag Kg/thm = ' 
yz27- Bosh Slag Kg/thm = ' 
yz28=' Lol in Lime (%) ' 

yz29='Heat Losses(cal/mole)= ' 
yz30~ 'Shaft F,fticiency(%) = ' 
yz3 1“' RAFT (K) 
yz32- Heat Demand(Mcal/thm)' 
yz33~'BIast Temperature(C)~ ' 
yz34~'Oxygcn enrichment (%)"= ' 
yz3 5- Reserve Zone Temp.(K)= ' 
yz36='Heat Losses(Gcal/thm)= ’ 
yz37~'Heat Losses(Pct ) =' 



^ ( 



c grO gangue to iron ratio in burden 
c snfo % of fe2o3 in sinter 
c rmnp . % recovery of manganese 
c fcp : % of carbon in hotmetal 
c fsip ‘ % of silicon in hotmetal 
c fmnp . % of manganese in hotmetal 
c fsp ' % of sulphur in hotmetal 
c fpp % of phosphorus in hotmetal 
c Hrhrn : heat of dissociation of silica, mno and p2o5 

c 

gor=por*(oi( 1 )+or(2)+or(3)+or(4)) 
gsr=( 1 00 0-por) * (sr( 1 )+sr(2)+sr(3 )+sr(4)) 
grO=(gor+gsry(por*or(7)+(1 00 0-por)*sr(5)) 
srfo=(sr(5)-sr(6)*0 78)* 1 43 
rmn=tmnp/ 100.0 
fc —fcp/ 1 00 0 
fsi— fsip/ 1 00 0 
fmn-fmnp/100 0 
fs - fsp/ 1 00.0 
fp =fpp/I00.0 
ffe= 1 -fc-fsi-fmn-fp-fs 
Wcl(I)=0.0 
Wcl(2)=2 0 
Wcl(3)-2.0 
Wcl(4)-2.0 
Wcl(5)=2.0 
c Wcl(I)~1.0 

c Wcl(2)=1.25 

c Wcl(3)=1.5 

c Wcl(4)=1.75 

c Wcl(5)=2.0 

rfl(I) =2073.0 
rft(2) =2073.0 
rft(3) "2073.0 
rft(4) =2073.0 
rft(5) =2073.0 
bt(l) =900 0 
bt(2) =900 0 
bt(3) =900 0 
bt(4) =900 0 
bt(5) - -900.0 
oxy(l) =2.0 
oxy(2) =2 0 
oxy(3) =2.0 



oxy(4) =2.0 
oxy(5) =2.0 
HIoss(l)-HIs 
Hloss(2)=HIs 
FIIoss(3)-nis 
Hloss(4)=Hls 
Hloss(5)=Hls 
efcy(l) =0.65 
efcy(2) =0 65 
efcy(3) =0 65 
efcy(4) =0 65 
efcy(5) ~0.65 
rzt(l) =1123 0 
rzt(2) =1123.0 
rzt(3) =1123 0 
rzt(4) =1 123 0 
rzt(5) =1123.0 
Mh2o(l) =0.005 
Mh2o(2) -0.005 
Mh2o(3) "0.01 
Mh2o(4) -0 015 
Mh2o(5) =0 02 
lmpt( I ) =0 0 
lmpt(2) =0 0 
lrnpt(3) ~ 1 0.0 
Impt(4) =15.0 
lmpt(5) =20.0 
i3=0 
f=0.5 
ack=0.01 
acl=0 0 1 
Loy =20.0 

do 70 j - 1,5 

lmpt(j)=0.00 

Mh2o(j)=0.01 

W1 =lmpt(j)*Wcl(j)/100.0 

His =Hloss(j) 

eta =efcy(j) 

Tbosh =rzt(j) 
tb-ht(j) t273,0 
tc= 1773.0 

Tbosha =Tbosh-298.0 
Tbsd =Tbosh-50.0 
Tbsda =I'bosha-50.0 



Lcaco3=(Loy/44 0) 

Lca=(ca*( 1-Loy/I00)) 

Lsi--(si*(l-Loy/I00)) 

Lmg=(mg*( 1 -Loy/ 1 00)) 

Lai- (al*( I -Loy/ 1 00)) 

o2cl Acl* Acl('c' t ' 1 5/1. 0c6 t ocl/3.2e3 

Wlcac=Lcaco3* 1 Oe-2 

WIsl= 1 .0-L,caco3 *0 44 

Wclsl-Acl*( 1 0-Aclfe*0 0l)*0.0I 

Wcksl Ack I 0-Ackle*0 01)*0.01 

Hrhm=fsi*HfSi/28 0+fmn*Hfmn/55 0+fp*Hfp/62.0 

Hrvvl Hfcac*Wlcac 

Hrwcl=HfFo3*Acl*Aclfe/I 6e6 +Hcrack(Ccl,Hcl,ScI,Ocl) 
c VARIABLES 

Wslcl=0 0 1 * Acl*( 1 .0-Aclfe*0 01) 

Wslck=0.0l *Ack*(l 0-Ackfe*0 01) 

WsIIi- I 0-Lcaco3*0 44 
oea=Oxy(j)*0 01 
02 =0.2 H oca 
N2 "1.0-02 

£ jfc 9fc ifc jfc Jfc ift jfc ♦ 

if (i3.ne. I ) go to 777 
779 continue 

Rcg= 0 001 12 

c variables 1 CO 2 SIO 3 SI02 4 Fe 5 COKE 
c CARBON BALANCE 
v( I , I )— 12.0 
v(l,2)=0.0 
v(l,3)=0 0 
v(I,4)=0.0 

v(I,5)= - (0.0l*Fcck) 
vr( 1 )=0.0 1 *Wcl(j)*Ccl+Wl*Lcaco3 *0.12 
c SILICON BALANCE 
v(2, 0=0.0 
v(2,2)" 1.0 
v(2,3)= 1.0 
v(2,4) - 0.0 

v(2,5)= - (Ack*Acksi/6 0E5) 
vr(2)=(Wcl(j)*AcI*AcIsi/6.0E5)fLsi*Wl/60 0 
c SiO BALANCE 
v(3,l)= 0.0 
v(3,2)= 1.0 
v(3,3)= 0 0 
v(3,4)= 0.0 

v(3,5)= •(f l 'Ack*Acksi)/6 0E5 





vi (3) ( f* Wcl(j ) * Acl * AcIsi)/0 UIi5 
c IRON BALANCE 
v(4, 1 ) 0 0 
v(4,?) 0 0 
vH.i) oo 
v(4,4) I 0 

v(4,5) -(2*Ack + Ackle)/l 6E6 
vi(4)-(2 t Wcl(j)*Acl*AcIfe)/1.6E6 
c OXYGEN BALANCE 
v(5, 1) - 1 0 
v(5,2)~-l 0 
v(5,3 ) 0 0 
v(5,4)= 0 0 

v(5,5)~ -(3*Ack*Ackfe)/ 1 6E6 

vi(5)- 2*02 1 (Wd(j)*0cl)/1 .6E3+0.02* Wl*Lcaco3+ 

* (3 * Wcl(j ) * Acl* Aclfe)/ 1 . 6E6+mh2o(j) 

Call SOLN(5,v,vr,s) 

c HYDROGEN BALANCE 

nh2-mh2o(j>(WcI(j)*Hd)/2.0E2 
c NITROGEN BALNCE 

nn2=N2 KWcl(j)+Wl)*Rcg+(Ncl*Wcl(j))/2.8E3 
c CaO BALANCE 

ncao=(Wl*Lca)/56.0 
c MgO BALANCE 

nmgo“(Wl*Lmg)/40.0 
c AI203 BALANCE 

nAl2o3~(Wd(j)*Acl*Aclal)/l.02E6+(s(5)*Ack*Ackal)/l.02E6+ 

* (Wi*Lal)/l 02E2 

c INPUT SENSIBLE HEATS 
c 1 1 EAT OF FORMATION OF CO GAS 

SreacI^sO )*Hfco 

c HEAT OF FORMATION OF SIO GAS 
Sreac2~-s(2)*(I ffsio) 
c SENSIBLE HEAT OF OXYGEN 
So2~02*cps(ac8,bc8,cc8,tb) 

c sensible; heat of nitrogen 

Sn2~N2*cps(ac6,bc6.,cc6,tb) 
c SENSIBLE HEAT OF WATER 
Sh2o=mh2o(j)*cps(ac9,bc9,cc9,tb) 
c SENSIBLE HEAT OF COKE 
Washck 0.0l*(S(5)*Ack) 

Scoke Washck*0.25*(tc-298)i-((S(5)*Fcck)/1.2E3)*cps(acI,bcI,ccl,tc) 

c OUTPUT SENSIBLE HEATS 

c HEAT OF FORMATION OF H20 

Sreac3~mh2o(j)*Hfh2o 



c SENSIBLE HEAT OF IRON 

Sfe- s(4)*(10 6*(tf- 1873) + 280 0*56 0) 
c SENSIBLE HEAT OF (Sio2+Al2o3+Cao+Mgo) 

Sslagl=(60 0*S(3)+102 0*nAl2o3+56.0*ncao+40 0*nmgo) 
c Sslag=SsIagl*(0 25*(tf-l873)+480) 
c I IEAT OF DISSOCIATION OF Fe2o3 
SFe2o3=(S(4)/2.0)*HfFo3 
c HEAT OF CRACKING 

Scrack-I Icrack(C.cl, I Icl, Scl.Ocl) * Wcl(j)+ W1 *Hrwl 
SLH=SreacI+Sreac2+So2+Sn2+Sh2o+Scoke-Sreac3+ 

* S(4)*(10 6*1575 0-280 0*56)- SslagI*(Hslag-0.25* 1575)- 

* SFc2o3-Sciack 

tfa=2000 
999 continue 
tf ~tfa 

Sh2 I= r nli2*cpa(ac4,bc4,cc4,tf) 

Sco l=s( 1 )*cpa(ac3,bc3,cc3,tf) 

Sn2I I _ nn2*cpa(ac6,bc6,cc6 1 tl) 

SsioL”s(2)*cpa(ac2,bc2,cc2,tf) 

SRH-SH21 1 Sco H-Sn2 1 H Ssio I+Sslag l *0 25+S(4)* 1 0.6 
tfa--(SLH/SRH)i298 
if (abs(tfa-tf) gt I 0) go to 999 
write(6,*) tfa 
rfl(j) “tf'a 
p(2) s( 5 ) 
goto 778 
777 continue 
tf~rft(j) 

c raceway variables I=CO 2= coke 3= H20 4= Fe 5= slag 6= SiO 
c CARBON BALANCE 
u(I,I)=I2.0 
u( l ,2)— -(0.01 *FCck) 
u(I,3)=0.0 
u( 1 ,4) "0.0 
u(I,5) 0.0 
u( 1,6) 0.0 

ur( 1 ) =0.0 1 * Wcl(j)*Ccl+ 12.0* Wl*Wlcac 
c OXYGEN BALANCE 
u(2,I)-1.0 

u(2,2)=-( Ack* Ackfe*3 0/( I 6E6)) 

u(2,3) -1.0 

u(2,4)-0.0 

u(2,5)=0.0 

u(2,6)~-I.O 

d I =WcI(j )*0.0 1 *(Ocl/ 1 6+(AcI* Aclfe*3 .0/1.6E4)) 



d2 (Wleac*Wh()2)*2 0 
ur(2) =dl ^(12 
SILICON BALANCE 
u(3,l)=0.0 

u(3, 2)~-(Ack*Acksi*f76.0E5) 

u(3,3)"0.0 

u(3,4)=0,0 

u(3,5)~0 0 

u(3,6)=I 0 

ur(3) =(Wcl(j)*Acl*Aclsi*t76 0E5) 

IRON BALANCE 
u(4, l)”0 0 

u(4,2)=-(Ack*Ackfe* 1 12 0/I.6E6) 

u(4,3)=0 0 

u(4,4)=l .0 

u(4,5)~0.0 

u(4,6)=0 0 

ur(4) -(Wcl(j)*Acl*Aclfe* 1.12/1 6E4) 

SLAG BALANCE 
u(5, 1 )“().() 
u(5,2) --Wcksl 
u(5,3)=0.0 
u(5,4)=0.0 
u(5,5)=I 0 
u(5,6)~60.0 

ur(5) ~WcI(j)*WclsI+Wl*WIsI 

HEAT BALANCE 

cpc =cps(acl,bcI,ccI,Tc) 

cpsio=cps(ac2,bc2,cc2,Tf) 

cpco =cps(ac3,bc3,cc3,Tf) 

cph2=cps(ac4,bc4,cc4,Tf) 

cpco2-=cps(ac5,bc5,cc5,Tf) 

cpn2“cps(ac6,bc6 ? cc6,Tf) 

cpo2tb~cps(ac8,bc8,cc8,Tb) 

cpn2tb-=cps(ac6,bc6,cc6,Tb) 

cph2o ~cps(ac9,bc9,cc9,Tb) 

cpfe ~0.2E-3*(Tf**2)+(Tf*9.77)-4636 79 

u(6,l)~-Hfco fcpco 

ur4 -Ack*2.5e~3*(Tc-298.0) 

u(6,2) ~-(ur4-t (FCck*cpc/1.2E3)) 

u(6,3)-Hfli2o+cph2-cph2o 

u(6,4) (Iinb3*0 5)/56.0t-280.0K(I0 56*(Tf-l873))/56.0) 
u(6,5) "Hslag i ( Tf-1873 0)/4.0 
u((),(i) I Ifsio i cpsio 

url _ -(WcI(j)’ ,, Hcrack(CcI,Hcl,Sd,Ocl)^WI*Hrwl) 
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ur2 =O2*cpo2tb+N2*cpn2tb-Wcl(j)*Hcl*cph2/200 0 
Mn2 ~N2 t (Wcl(j)*Ncl/2800 0)4((Wci(j)+Wl)* 1.12E-3) 
ur(6) =uri+ur2-Mn2*cpn2 
call S0LN(6,u,ur,p) 
if(p(3).It 0 0) then 
Mh2o(j)- 0 0 
go to 779 
else 

Mh2o(j)~p(3) 

endif 

778 continue 


c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SECOND S TAGE MODELLING 


Mbh2 ■ moles of hydrogen 
Mbn2 moles of nitrogen 
cpo2tb ' specific heat of oxygen at Tb temperature 
cpn2tb specific heat of nitrogen at Tb temperature 
Wsit : weight of silica at tuyere/ton of hotmetal 
p2 . weight of coke 
y(5) : weight of hotmetal 
cat : weight of calcium at tuyere 
alt : weight of alumina at tuyere 
FOin . moles of iron as fe2o3 
FEOin . moles of iron as feo 
Wbd : weight of burden /ton of hotmetal 
Wor weight of ore /ton of hotmetal 
Wsn ' weight of sinter /ton of hotmetal 
Wcoke . weight of coke /ton of hotmetal 
Wcoal • weight of coal /ton of hotmetal 
Wftiel weight of total fuel 
Wslib . weight of slag /'ton of hotmetal 
BSsi weight of silica in bosh 
BSca ' weight of calcium in bosh slag 
BSmg . weight of magnesium in bosh slag 
BSal : weight of al2o3 in bosh slag 
Wbs . weight of bosh slag 
bs3 : basicity of slag in stack 
bs4 : basicity of slag in stack 
cka : weight of ash in coke 
cla . weight of ash in coal 
CSsi . weight of silica in cast slag 
CSca weight of calcium in cast slag 
CSmg ■ weight of magnesium in cast slag 
CSal . weight of alumina in cast slag 
CSmn weight of manganese in cast slag 





c 

sib 

% of silica in bosh slag 

c 

cab 

% of calcium in cast slag 

c 

mgb 

% of magnesium in cast slag 

c 

bs5 

basicity in bosh slag 

c 

bs6 

basicity in bosh slag 

c 

slsi 

. % silica in cast slag 

c 

slca 

: % calcium in cast slag 

c 

slmg 

% magnesium in cast slag 

c 

slal 

% aluminium in cast slag 

c 

sis 

% sulphur in cast slag 

c 

slmn 

% manganese in cast slag 

c 

slfe 

% feo in cast slag 

c 

bsl 

basicity in cast slag 

c 

bs2 

• basicity in cast slag 


c variable l=Coke 2 = FeO 3= Fe 4= CO 5= hot metal 6= slag 
c Oxygen balance in lower zone of BF 
a(l, !)=-(! 5*Ack*Ackfe)/I 6E6 
a(l,2)=-0 5 
a( 1 ,3 )"0 0 
a( 1 ,4)=0 5 

a(I,5)=-(fsi/28 +fmn/l l0.+fp*2.5/62 ) 
a(l,6)=ffeos*0 01/144. 
c Iron Balance in lower zone of BF 
a(2, 1 )=-Ack* Ackfe*0 7E-4 
a(2,2)=-56 0*0 95 
a(2,3)=-56.0 
a(2,4)-=0 0 
a(2,5)=ffe 

a(2,6)=0 95*0.56*ffeos/72 0 
c Carbon balance in the lower part of BF 
a(3,l)=-FCck/I00 
a(3,2)=0.0 
a(3,3)=0.0 
a(3,4)=I2. 
a(3,5)=fc 
a(3,6)-0 0 
c Slag Balance 
a(4,I)=-Wcksl 
a(4,2)=-gr0*56.0*0.95 
a(4,3)~-gr0*56.0 
a(4,4)=0 0 

a(4,5)=(2 14*fsi-(I-rmn)* 1 29*fmn/rmn) 
a(4,6)=I 0-ffeos*0.01-fs*ks 
c Shaft Equilibrium FeO + CO = Fe + C02 





a(5,l)=0 0 
a(5,2)=0 0 
a(5,3)= 1/0.95 

K2=10**(-(4350 0/(4 606*Tbsd))+0 5 1) 

K.l=I0**((4300.0/(4.606*Tbsd))-l 14) 

a(5,4)=-eta*Kl/(l+Kl) 

a(5,5)-0.0 

a(5,6)=0 0 

c Heat balance in the Lower Zone of BF 

a(6, 1 )~-(cps(ac 1 ,bc 1 ,cc I ,'I bsd)*FCck/ 1 2.0+Ack*0 25*Tbsda)*0 0 1 
1 +Ack*Ackfe*hffe/I 6e6 

a(6,2)“-cps(ac7,bc7,cc7,Tbsd)t HfTeo-0 25*Tbsda*gr0*56.0*0 95 
a(6,3 )=-ac9 *Tbsda-0 25*Tbsda*56 0 + gr0 
a(6,4)=cp s(ac3 ,bc3 , cc3 ,Tbo sh)-Hfco 

a(6,5)=Hhm+Hrhm-0 25*(142*fp/62+(I .-rmn)*I 29*fmn/rmn+fs)*Tbsda 
a(6,6)~HsIag~ffeos*HfFeo/7200 0-0 25*Tbsda*fs*ks 
SHh2=cps(ac4,bc4,cc4,Tbosh) 

SHn2=cps(ac6,bc6,cc6,Tbosh) 

cpc =cps(acI,bcI,ccl,Tc) 

cpsio=cps(ac2,bc2,cc2,Tf) 

cpco =cps(ac3,bc3,cc3,Tf) 

cph2=cps(ac4,bc4,cc4,Tf) 

cpco2-cps(ac5,bc5,cc5,Tf) 

cpn2=cps(ac6,bc6,cc6,Tf) 

cpo2tb=cps(ac8,bc8,cc8,Tb) 

cpn2tb=cps(ac6,bc6,cc6,Tb) 

cph2o =cps(ac9,bc9,cc9,Tb) 

cpfe =0 2E-3*(Tf**2)+(Tf*9.77)-4636.79 

Mn2 ~N2H (Wcl(j)*Ncl/2800.0)-t ((Wcl(j)+WI)* 1 12E-3) 

Mbh2=(Wcl(j)*Hcl)/200 0+Mh2o(j) 

Mbn2=((Wcl(j)+Wl)* 1, 12E-3)+(Wcl0*Ncl/28OO 0)+0.79-0ea 
c THE COEFFICIENTS OF THE MATRIX 
c( I )= W1 * W lcac+ W cl(j) * o2cl+02+mh2o(j)/2 . 0 
c(3)=Wcl(j)*Ccl*0 0I+Wl*Wlcac*I2 
c(2)~Wcl(j) + Acl*Aclfe*0 7e-4 
c(4)=WcI(j)*Wclsl 
c(5)=Mbh2*eta*K2/(I 0+K2) 
c6a=02*cpo2tb 
c6b~N2*cpn2tb 
c6c=Mh2o(j)*(cph2o-Hfh2o) 

c6d=Mbh2*SHh2+Mbn2*SHn2+Hls*0 8+WlsI*WI*0.25*Tbsda 
c(6)~c6a tc6b+c6c-WcI(j)*HrwcI-Wl*HrwI-c6d 
call SOLN(6,a,c,y) 

call bdn(por,or(7),y(2),y(3 ),sr(5 ),srfo,sr(6),h) 



SLAG FORMED IN THE TUYERE 
Wsit=(p(2)*WcksI+Wcl0*WcIsI+Wl*WlsI)*lOOO 0/y(5) 
sit =^p(2)*Ack*Acksi+WcI(j)*AcI*AcIsi+Wl*Lsi* 10000 0)*0 I/y(5) 
cat =Wl*Lca*I000 0/y(5) 
mgt =WI*Lmg*I000 0/y(5) 

alt =((p(2)*Ack*AckaI+WcI(j)*AcI*AcIaI+WI*LaI* 10000 0)*0 I)/y(5) 
write(6,*)alt 
FOin=h(4)* I000.0/y(5) 

FEOin=h(5)*IOOO 0/y(5) 

Wbd =h(I)*1000 0/y(5) 

Wor =h(6)*1000 0/y(5) 

Wsr =h(7)*I000.0/y(5) 

Wcoke(j)=y( I ) * 1 000 0/y(5) 

WcoaI0=WcI0* 1 000 0/y(5) 

WfueI0=Wcoke(j)+Wcoai(j) 

WIime0=WI* 1 000 0/y(5) 

WsIag0=IOOO 0*y(6)/y(5) 

W sIib(j)=W sIag(j)-W sit 
siuib=60000*fsi/28.0 
mnuib=7 1 000*fmn/ 5 5 0 
BIast0=224OO. 0*( I O+Mh2o0)/y(5) 

SLAG COMING TO THE BOSH FROM THE STACK 

B S si=( W or * or( 1 )+W sr* sr( I )) * 0 01-Wlime(j)*Lsi 

B S ca=(W or* or(2)+W sr * sr(2 )) *0. 0 1 - Wlime(j) *Lca 

BSmg=^or*or(3)+Wsr*sr(3))*0.0I-Wlime(j)*Lmg 

BSaI=(Wor*or(4)+Wsr*sr(4))*0 0I-WIime(j)*Lai 

Wbs(j)=BSsi+BSca+BSmg+BSaI 

bsi(j)=BSsi* 1 00. 0/Wbs(j) 

bca0=BSca* IOO.OAVbs0 

bmg(j)=BSmg*lOO.O/Wbs0 

baI(j)=BSal* 1 00. 0/Wbs(j) 

bs3(j)=BSca/BSsi 

bs4(j)=(BSca+BSmg)/BSsi 

cka =W coke(j ) * Ack* 0 . 0 1 

cla =WcoaI0*AcI*O.Ol 

CSsi=(Wor*or(I)+Wsr*sr(l)+cka*Acksi+cIa*Acisi)*0.0I-fsi*2.IE3 
CSca={W or*or(2)+Wsr*sr(2))*0. 0 1 
CSmg=(Wor*or(3)+Wsr 3 * : sr(3))*O.OI 

CSaI^Wor*or(4)+Wsr*sr(4)+cka*AckaH-cIa*AcIaI)*0.0I 

CSmn=( 1 . 0-rmn)*fmn* 1 290. 0/rarn 

Wcs(j)=(CSsi+CSca+CSmg+CSaI+CSmn)/(I 0-ffeos*0.01-fs*ks) 
SLAG FORMING AT THE BOSH 
sib(j) =(CSsi-sit+siuib)* I0O.0AV siib0 
cab0 =(CSca-cat)*lOO.O/WsIib0 
mgb(j) =(CSmg-mgt) * 1 00. OAV slib(j) 



alb(j) =(CSal-alt)* 100 0/WsIib(j) 
write(6,*)alb 

febO') =fFeos*WsIag(j)AVsIib(j) 
sb 0 =ks*fsp*WsIag(j)AVsIib(j) 
mnb(j) =(CSmn+mnuib)* 100 O/Wslib0 
bs5(j) =cab(j)/sib(j) 
bs6(j) =(cab0+mgb0)/sib(j) 

^ *************** 

sIsi(j)=C Ssi * 1 00 . 0/W cs(j) 
sIca(j)=CSca*100 OAVcs(j) 
stmg(j>CSmg*IOO OAVcs(j) 
sIaI0=CSaI* 100 OAVcs(j) 
sls(j) =ks*fsp 

slmn0=CSmn* 1 00. 0/WcsQ 
sIfe(j)=ffeos 
bsl(j) -slca(j)/slsi(j) 
bs2(j) =(slca0+slmg(j))/ slsi(j) 
ph2o(j)=Mh2o(j)* 1 .8E4/22 4 
c P h2o(j)=Mh2o0* 1 ,8E4/((I ,O+Mh2o0)*22.4) 

c THIRD STAGE MODELLING 

c 

c CL • weight of coal /ton of hotmetal 
c CK : weight of coke /ton of hotmetal 

c LM : weight of lime /ton of hotmetal 

c SL : weight of slag /ton of hotmetal 
c xms : % moisture in coke on dry basis 

c msck : moisture in coke in moles 

c 02ty : moles of oxygen /ton of hotmetal 

c N2ty : moles of nitrogen /ton of hotmetal 

c ffes : % of iron in slag 
c Mst ‘ moles of H2o /ton of hotmetal 
c feci : moles of iron in coal 
c feck : moles of iron in coke 
c FEot : moles of iron out 
c FEin : mloes of iron in 
c CALm : moles of cacao3 
c Tco2 : moles of co2 
c Tco : moles of co 
c Th2o : moles of water 
c Th2 : moles of hydrogen 
c t : total moles of gases 

c 

CL=WcoaI0 

CK=Wcoke0 

LM=WIime(j) 



SL=WsIag(j) 

BL=1000 0/y(5) 

msck=0 0l*xms*CK/(I8 0*(I 0-0 0l*xms)) 

02ty=BL*02 
N2ty=BL*N2 
fFes=0 56*fFeos/72 0 
Mst=BL*Mh2o(j) 

Tn2=Mn2*BL 
FEci=CL* Acl* AcIfe/0. 8E6 
FEck=CK* Ack* Ackfe/0 . 8E6 
FEot=(ffe *1000. O+fFes* SL)/5 6 0 
FEin==FEot-FEcI-FEck 
CAim=LM*Lcaco3 * 1 . 0E-4 

Cin(CL*Ccl+CK*FCck)*0 01+I2.0*CAlm-fc*I000 0)/I2 0 
C12=(CL*OcI/I 6E3)+2 0*BL*02+Mst+2.0*CAIm+(FEcl+FEck)*1.5 
C22=h(3)*FEin+(fsi/14 0+fmn/55 0+fp/1 2.4)* 1000 0 
C2=CI2+C22 
C3=Mst+(CL*HcI/200. 0) 

C4=EXP((8650 0/773 . 0-7. 59)/ 1 . 9) 

C5C4-1.0 

C6=C2-C4*C2-C I-C3 
C7=C4*C I *(C2-C I ) 

Tco2=root(C5,C6,C7,Cl) 

Tco =Cl-Tco2 

Th2o=C2-C I -Tco2 

Th2 =C 1 +C3 +Tco2-C2 

t=Tco+Tco2+Th2+Tn2+Th2o 

C8 =(t-Th2o)*0 01 

tt(j)=t*22.4 

qa(j)=Tco/C8 

qbO>Tco2/C8 

qc0=Th2/C8 

qdO'HTco2* 1 00. 0)/(T co+T co2) 
qe(j)=Tn2/C8 

Hin=Hfco*Tco+Hfco2*Tco2+02ty*cpo2tb+N2ty*cpn2tb+Mst*cph2o 
HI =Hfh2o * (Mst-Th2o )+Hfcac*C Alm+SL *HsIag+ 1 000. 0 *Hhm 
H2 ={Hfsi*fsi/28 . 0+Hfmn*fmn/ 5 5 0+Hfp*fp/62.O)*I00O.0 
H5 ={FOin+0. 5 *(FEcI+FEck))*Hffo3+FEOin*Hffeo 
H3 =H5-(ffes*SL*HfFeo/56.0) 

H4 =Hcrack(CcI,HcI,Scl,Ocl)*CL+HIs*BL+msck*(Hvp+cph2o) 

Stpg=Hin-H I -H2-H3 -H4 

Tg=0.0 

Tgp=700.0 

itg=I 

100 DF(abs(Tg-Tgp).Ie.3.0) GO TO 101 



1 0 s 


if(itg ge 20)go to 10 1 
itg=itg+I 

Tg=Tgp 

cco =Tco*cf(ac3,bc3,cc3,Tg) 
cco2=Tco2*cf(ac5,bc5,cc5,Tg) 
ch2 =Th2*cf(ac4,bc4,cc4,Tg) 
ch2o=Th2o * cf(ac9,bc9,cc9 ,Tg) 
cn2 =Tn2*cf(ac6,bc6,cc6,Tg) 

Tgp =Stpg/(cco+cco2+ch2+ch2o+cn2)+298 0 
GO TO 100 
101 Tgt(j)=Tg-273.0 

Hdl =Tco*cps(ac3,bc3,cc3,Tg)-KTco2*cps(ac5,bc5,cc5,Tg) 
Hd2 =Th2 * cp s(ac4 ,b c4, cc4,T g)+Tb2 o * cp s( ac9,bc9,cc9,Tg) 
Hd3 =Tn2*cps(ac6,bc6,cc6,Tg) 

Hdem(j)= 1 . 0E-3 *(H1+H2+H3+H4+Hdl+Hd2+Hd3) 

tI=mh2o(j) 

tk=tb-273 

call RAF(tk,02,tI,L) 
write(6,*) L 
tll=ph2o(j) 

to= 60 . 0 *( 0.0 1 *oxy(j)/(0- 7 9-0.01 *oxy(j))) 
tpc=((WcI(j)*60 0)/ (22 4*1.0E3)) 
call RAFT (tk,tl 1 ,to,tpc,L 1 ) 
write(6,*) LI 

Gtm(j)=(HIs*BIast(i)*I0E-6)/22.4 

PGtm0=((HIs*BIast0)/(HdemQ*224 0)) 

70 CONTINUE 
rr(I)=0.0 

do 71 j=2,jt 

rr0=(Wcoke( I )-Wcoke0)/W coal(j) 

71 continue 
write(21,80) 

80 fonnat(20x,'COMPirTED OPERATING VARIABLES') 
write(21,8I) 

81 format(I9x,' >!f ) 

writ e(2 1,39) 

39 format(lx,80(Ih_)//) 
c write(21,98)yz28 

98 format(45x,a20) 
c write(2 1, 83 )(Loi(j) j == I Jt) 

write(21,84)yz3 I,(rft(j) j®* J*) 
c write(2I,201)yz31 
c write(21,83)(rft0j=IJt) 

write(21,84)yz33,(bt(j)J=I,jt) 
c write(21,20I)yz33 



! 0 


c write(2 1 ,83)(bt(j),j= I jt) 

write(2 1 ,84)yz34,(oxyQj= 1 jt) 
c write(2I,20I)yz34 
c write(2 1 ,208)(oxy(j),j= I jt) 
write(2 1 ,84)yz3 5,(rzt(j) j= I ,jt) 
c write(2I,20I)yz35 
c writ e(2 1,83 )(rzt(j),j= l,jt) 

208 format(/30x,5(f7.2,3x)/) 

write(2 1 ,84)yz29,(HIoss(j) j= I Jt) 
c write(2I,20I)yz29 
201 format(45x,a22) 
c write(2I,83)(HlossQj=Ijt) 

write(2 1 ,84)yz3 0,(efcy(j)* 1 00 ,j= I jt) 
c write(2I,20I)yz30 
c write(2 1 ,83 )(efcy(j)* 1 00 j= I jt) 

write(2 1 ,84)yz2,(WcoaI(j) j= I jt) 
writ e(2 1,39) 
c write(21,82)yz2 
c82 format(45x,a20) 
c write(21,83)(WcoaI(j)j=I jt) 
c83 format(/30x,5(f7. 1,3x)/) 

write(2I,84)yyI,(ph2o(j)j=I jt) 

84 format(/5x,a23,5x,5(f7.2,3x)) 
write(2 1 ,84)yz I ,(Wcoke(j) j= I jt) 
write(2 1 , 8 5)yz 1 3 ,(rr(j) j= 1 jt) 
write(2 1 ,84)yz6,(WIime(j) j= I jt) 

85 format(/5x,a20,5x,5(f7 2,3x)) 
write(2 1 ,84)yz3 ,(WfueI(j),j= I jt) 
write(2 1 ,84)yz4,(W slag(j) j=I ,jt) 
write(2 1 ,84)yz25,(Wbs(j) j= I jt) 
write(2 1 ,84)yz27,(WsIib(j)j= I jt) 
write(2 1 ,84)yz5,(BIast(j) j= 1 jt) 
write(2 1 ,84)yz3 6,(Gtm(j) j= I jt) 
write{2 1 ,84)yz37,(PGtm(i) j = I jt) 
write(2 1,84)yz I2,(tt(j) j=I jt) 
write(21,84)yzl I,(Tgt(]),j=I jt) 
write(2 1 ,86)yz32,(Hdem(j) j= I jt) 

86 format(/5x,a22,4x,5(f7. 1,3x)) 
write(2I,9I) 

91 format(/5x/CAST SLAG COMPOSITION ) 
writ e(2 1,92) 

92 format(5x/ ) 

write(2 1 ,85)yz 1 6,(slsi(j),j= I jt) 
write(2 1 ,85)yz 1 7,(slca(j) ,j= 1 jt) 
write(2 1 ,85)yz 1 8,(slmg(j) ,j= 1 jt) 



write(2 1 ,85)yz 1 9,(sIaI0 j= I ,jt) 
write(2 1 ,85)yz24,(sis(j) j= Ljt) 
write(2 1 ,85)yz2O,(sIfe0 j= I ,jt) 
write(2 1,85)yz2 1,(simn(j) j= I jt) 
write(2 1 ,85)yz22,(bs I 0 j= I jt) 
write(2 1 ,85)yz23 ,(bs20 j= I ,jt) 
write(2I,I3) 

1 3 format(/5x,'STACK SLAG BASICITY *) 
write(2I,I4) 

14 format(5x, r *) 

write(2 1 ,8 5)yz 1 6, (bsi(j),j= 1 jt) 
write(2 1 ,85)yz 1 7,(bca(j) j= I jt) 
write(2I,85)yzI8,(bmg0j=ijt) 
write(2 1 ,85)yz 1 9, (bal(j ) j= I jt) 
write(2 1 ,8 5)yz22, (bs3 (j) j= I jt) 
write{2 1 ,85)yz23 ,(bs40 j= I jt) 
write(2I,I5) 

1 5 format(/5x,BOSH SLAG BASICITY j 
write(21,2I) 

21 format(5x/ •*) 

write(2 1,85)yzI6,(sib(j)j=I,jt) 
write(2 1 ,85)yz 1 7 ,(cab(j) j= I jt) 
write(2 1 ,85)yz 1 8,(mgb(j)j= I jt) 
write(2 1 ,85)yzi9,(alb0 j= I jt) 
writef 2 1 ,85)yz24,(sb0 j= I jt) 
write(2I,85)yz20,(feb(j)jf= I jt) 
write(2 1 ,85)yz2 l,(mnb0, j= I jt) 
write{2I,85)yz22,(bs50,j=T jt) 
write(21,85)yz23,(bs60j=Ijt) 
writ e(2 1,89) 

89 format(/5x,TOP GAS COMPOSITIONCdry basis)') 
writ e(2 1,90) 

90 format(5x/ j 

write(2 1 ,8 5)yz7 ,(qa(j) ,j= 1 jt) 
write(2I,85)yz8,(qb(j)j = I jt) 

write(2 1,85)yzI4,(qd0,j=Ijt) 
write<2 1,85)yz9,(qc(i)j= I jt) 
write(21,85)yzl 0,(qe(j) j=l jt) 
write(2I,72) 

72 format(5x,40(lH*)//) 

410 continue 
stop 
end 

SUBROUTINE SOLN(me,u,c,y) 
dimension u(10,I0),c(I0),y(I0) 



mej=me+ 1 
do I00im=I,me 
1 00 u(Im,mej)=c(im ) 
do 1 15 i=l,me 
do 1 05 j= I, meg 
if[j EQ i) go to 105 
u(i.j)=u(i,j)/'u(i,i) 

105 continue 
u(i,i)=I 

do 1 15 ii=I,me 
if(i EQ ii) go to 1 1 5 
do I10j=i+I,mej 

1 10 u(ii,j)=u(ii j)-u(i,j)*u(ii,i) 

115 continue 

do 120im=I,me 
120 y(im)=u(ira,mej) 
return 
end 

SUBROUTINE bdn(a,n,b,c,d,e,f,x) 
real g,xI2,x22,a,b,c,d,e,£n,x(10),wo,ws 
g={a*n+(I00.0-a)*d)/l ,0e4 
x(l)=56 0*(b+c)/g 
ws =x(l)*0.01*(100.0-a) 
wo =x(I)*0 0I*a 
xI2=e*3. 0/160 0+&72.0 
x22=e/80 0+E72.0 
x(2)=xI2/x22 

x(3 j=(d*ws*x(2)+n*wo* 1 . 5)/(d*ws+n*wo) 
x(4)=x(I)*((100.0-a)*e/I60.0+a*n/T12 0)*1 0e-4 
x(5)=x(I)*(I00 0-a)*f77.2e5 
x(6)=wo 
x(7)=ws 
return 
end 

SUBROUTINE RAF (a,b,c,G) 
real G,a,b,c 

G=((0.934I*a+-8208.0*b-c*(2402-I.2I77*a)+94 76)/(I 0+b+c)) 

return 

end 

SUBROUTINE RAFT(a,b,c,d,GI) 
real Gl,a,b,c,d 

GI=(I559.0+0 839*a-6 033*b+82.87*c-0 03I*(I6666*d)**2) 
return 
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REAL FUNCTION cps(a,b,c,t) 
real tl 
1 1=298 0 

cps=(M:l)*(aH>*(t+tI)*0 5-c/(t*tI)) 

return 

end 

REAL FUNCTION cf(a,b,c,t) 
real tl 
1 1=298 0 

cf=a+b*(t+tl)*0 5-c/(t*ti) 

return 

end 

REAL FUNCTION cpa(a,b,c,t) 
realt 

cpa=a+b*(t+298)*0 5-c/'(t*298) 

return 

end 

REAL FUNCTION root(a,b,c,e) 
real d,x,y 

cf=SQRT(b*b-4 0*a*c) 
x=(-b+d)/(2 0*a) 
y=(-b-d)/(2 0*a) 
if(x .LE e) THEN 
root=x 
else 
root=y 
endif 
return 
end 

REAL FUNCTION Hcrack(a,b,c,d) 
real Hfco2,Hfh2o,hfso2 
c Hfco2 =94260 0 

Hfh2o =58850.0 
c Hfso2 =218334.0 

c Hcrack =0.0 1 *(a*Hfco2/I2 0+b*Hfh2o/2.0+c*ffiso2/32.0)-d 

c Hcrack =0.01 *a*5533 0/12 0 

Hcrack = 0.0l*(d*Hfh2o/I6 0+(b-d/8 0-c/I6.0)*2766 5) 

return 

end 



APPENDIX: DATAFILE 


!c& 


17 0 10 0 30 0 60 0 5 3 I 62 3 9 71 8 0 38 
19 0 10 0 300 60 0 81 0 40 
0 0236 0 9467 0 02476 0 00484 4 0 0 001 

26420.0 32380 0 168980 0 117540 0 40250 0 7250 9 
4 03 1.04E-3 2.04E5 

7.7 0.74E-3 0 70E5 
6 79 0 98E-3 0.11E5 
6 52 0 78E-3 0.12E5 

10 57 2.10E-3 2.06E5 

6 83 0 90E-3 0.12E5 

11 66 2.00E-3 0.67E5 
7.16 1.00E-3 0.40E5 

7 3 2.46E-3 0.0 

57800.0 196800.0 217300 0 92050 0 367060 0 134510 0 
63500 0 480 0 280 0 26420 0 5533 0 

195450.0 193000.0 94260 0 30.0 

4 19 0 95 0 1 0217 0 061 70.0 20 19 
1100 2273 0 1073 0 2000 0 5 2.2 1 
33.3 0 65 3 75 

1.7 0.2 0.1 1.2 0 053 0 3 68 95 
4.1 10 04 2.36 2.59 56.52 9 28 

10.0 0 28 100 
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